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(Dated: February 5, 2008) 

Abstract 

We describe the status of our effort to realize a first neutrino factory and the progress made in 
understanding the problems associated with the collection and cooling of muons towards that end. 
We summarize the physics that can be done with neutrino factories as well as with intense cold 
beams of muons. The physics potential of muon colliders is reviewed, both as Higgs Factories and 
compact high energy lepton colliders. The status and timescale of our research and development 
effort is reviewed as well as the latest designs in cooling channels including the promise of ring 
coolers in achieving longitudinal and transverse cooling simultaneously. We detail the efforts being 
made to mount an international cooling experiment to demonstrate the ionization cooling of muons. 

PACS numbers: 13.10.+q, 14.60.Ef, 29.27.-a, 29.20.Dh 
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I. INTRODUCTION 

Recent results from the SNO collaboration Q| coupled with data from the SuperK col- 
laboration [3] have provided convincing evidence that neutrinos oscillate and that they very 
likely do so among the three known neutrino species. Experiments currently under way or 
planned in the near future will shed further light on the nature of these mixings among neu- 
trino species and the magnitudes of the mass differences between them. Neutrino oscillations 
and the implied non-zero masses and mixings represent the first experimental evidence of 
effects beyond the Standard Model, and as such are worthy of vigorous scientific study. 

This document indicates our progress along a path toward establishing an ongoing pro- 
gram of research in accelerator and experimental physics based on muon beams, and neutrino 
beams derived therefrom, that can proceed in an incremental fashion. At each step, new 
physics vistas open, leading eventually to a Neutrino Factory and possibly a Muon Collider. 
This concept has aroused significant interest throughout the world scientific community. In 
the U.S., a formal collaboration of some 110 scientists, the Neutrino Factory and Muon 
Collider Collaboration, also known as the Muon Collaboration (MC) Q], has undertaken 
the study of designing a Neutrino Factory, along with R&D activities in support of a Muon 
Collider design. The MC comprises three sponsoring national laboratories (BNL, FNAL, 
LBNL) along with groups from other U.S. national laboratories and universities and indi- 
vidual members from non-U. S. institutions. 

One of the first steps toward a Neutrino Factory is a proton driver that can be used 
to provide intense beams of conventional neutrinos in addition to providing the intense 
source of low energy muons (from pion decay) that must first be "cooled" before being ac- 
celerated and stored. Our vision is that while a proton driver is being constructed, R&D 
on collecting and cooling muons would continue. A source of intense cold muons could 
be immediately used for physics measurements, such as determining the electric and mag- 
netic dipole moments of the muon to higher precision, muonium-antimuonium oscillations, 
muon spin rotation experiments and rare muon decays. Once the capability of cooling and 
accelerating muons is fully developed, a storage ring for such muons would serve as the 
first Neutrino Factory. Its specific beam energy and its distance from the long-baseline 
experiment will be chosen using the knowledge of neutrino oscillation parameters gleaned 
from the present generation of solar and accelerator experiments (Homestake, Kamiokande, 
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Super Kamiokande, SAGE, GALLEX, K2K, SNO), the next generation experiments (Mini- 
BooNE, MINOS, CNGS, KamLAND, Borexino), and the high-intensity conventional beam 
experiments that would already have taken place. 

A Neutrino Factory provides both v ^ and V e beams of equal intensity from a stored \x~ 
beam, and their charge-conjugate beams for a stored /x + beam. Beams from a Neutrino 
Factory are intense compared with today's neutrino sources. In addition, they have smaller 
divergence than conventional neutrino beams of comparable energy. These properties permit 
the study of non-oscillation physics at near detectors, and the measurement of structure 
functions and associated parameters in non-oscillation physics, to unprecedented accuracy. 
Likewise, they permit long-baseline experiments that can determine oscillation parameters 
to unprecedented accuracy. 

Depending on the value of the parameter sin 2 2# 13 in the three-neutrino oscillation for- 
malism, the oscillation v e — > is expected to be measurable. By comparing the rates for 
this channel with its charge-conjugate channel V e — > V^, the sign of the leading mass dif- 
ference in neutrinos, 5m| 2 , can be determined by observing the passage through matter of 
the neutrinos in a long-baseline experiment. Such experiments can also shed light on the 
CP-violating phase, 5, in the lepton mixing matrix and enable the study of CP violation 
in the lepton sector. (It is known that CP violation in the quark sector is insufficient to 
explain the baryon asymmetry of the Universe; lepton sector CP violation possibly played 
a crucial role in creating this asymmetry during the initial phases of the Big Bang.) 

While the Neutrino Factory is being constructed, R&D aimed at making the Muon Col- 
lider a reality would be performed. The Muon Collider will require muon beams that are 
more intensely cooled and have generally more challenging properties than those for a Neu- 
trino Factory, so the latter forms a practical goal en route to the former. A Muon Collider, 
if realized, provides a tool to explore Higgs-like objects by direct s-channel fusion, much as 
LEP explored the Z. It also provides a potential means to reach higher energies (3-4 TeV 
in the center of mass) using relatively compact collider rings. 



A. History 

The concept of a Muon Collider was first proposed by Budker jj| and by Skrinsky in 
the 60s and early 70s. However, additional substance to the concept had to wait until the 



6 



idea of ionization cooling was developed by Skrinsky and Parkhomchuk 6 . The ionization 

o n 

cooling approach was expanded by Neuffer |7j and then by Palmer j8( , whose work led t o th e 
formation of the Neutrino Factory and Muon Collider Collaboration (MC) 0] in 1995 162 1. 

The concept of a neutrino source based on a pion storage ring was originally considered by 
Koshkarev (l2|. However, the intensity of the muons created within the ring from pion decay 
was too low to provide a useful neutrino source. The Muon Collider concept provided a way 
to produce a very intense muon source. The physics potential of neutrino beams produced 
by high-intensity muon storage rings was briefly investigated in 1994 by King |l3(and in 
more detail by Geer in 1997 at a Fermilab workshop [lJ . lis] ] where it became evident 
that the neutrino beams produced by muon storage rings needed for the Muon Collider 
were exciting in their own right. As a result, the MC realized that a Neutrino Factory 
could be an important first step toward a Muon Collider. With this in mind, the MC has 
shifted its primary emphasis toward the issues relevant to a Neutrino Factory. The Neutrino 
Factory concept quickly captured the imagination of the particle physics community, driven 
in large part by the exciting atmospheric neutrino deficit results from the SuperKamiokande 
experiment. The utility of non-oscillation neutrino physics from neutrinos produced by muon 
storage rings has been studied in detail from 1997 onwards ^(|. 

There is also considerable international activity on Neutrino Factories, with international 
conferences held at Lyon in 1999 Q, Monterey in 2000 0, Tsukuba in 2001 Q, Lon- 
don in 2002 2lJ and another planned in New York in 2003 j^. There are also efforts 
in Europe [22 1 and Japan 2^] to study different approaches to realizing the neutrino fac- 
tory. Recently a proposal has been submitted to perform an International Muon Ionization 
Cooling Experiment (MICE) to the Rutherford Appleton Laboratory 24 1 . 



B. Feasibility Studies 

Complementing the MC experimental and theoretical R&D program, which includes work 
on targetry, cooling, rf hardware (both normal conducting and superconducting), high-field 
solenoids, liquid hydrogen absorber design, muon scattering experiments, theory, simula- 
tions, parameter studies, and emittance exchange (2^, the Collaboration has participated 
in several paper studies of a complete Neutrino Factory design. 

In the fall of 1999, Fermilab, with help from the MC, undertook a Feasibility Study 
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("Study-I") of an entry-level Neutrino Factory Study-I showed that the evolution of the 
Fermilab accelerator complex into a Neutrino Factory was clearly possible. The performance 
reached in Study-I, characterized in terms of the number of 50-GeV muon decays aimed at 
a detector located 3000 km away from the muon storage ring, was N = 2 x 10 19 decays per 
"Snowmass year" (10 7 s) per MW of protons on target. 

Simultaneously, Fermilab launched a study of the physics that might be addressed by 
such a facility 27[ and, more recently, initiated a study to compare the physics reach of a 
Neutrino Factory with that of conventional neutrino beams j3| powered by a high-intensity 
proton driver (referred to as "superbeams"). As will be described later in this paper, a 
steady and diverse physics program will result from following the evolutionary path from a 
superbeam to a full-fledged Neutrino Factory. 

Subsequently, BNL organized a follow-on study ( "Study-II" ) on a high-performance 
Neutrino Factory again in collaboration with the MC. Study-II demonstrated that BNL was 
likewise a suitable site for a Neutrino Factory. Based on the improvements in Study-II, the 
number of 20-GeV muon decays aimed at a detector located 3000 km away from the muon 
storage ring, was N — 1.2 x 10 20 decays per Snowmass year per MW of protons on target. 
Thus, with an upgraded 4 MW proton driver, the muon decay intensity would increase to 
4.8 x 10 20 decays per Snowmass year. (R&D to develop a target capable of handling this 
beam power would be needed.) Though these numbers of neutrinos are potentially available 
for experiments, in the current storage-ring design the angular divergence at both ends of 
the production straight section is higher than desirable for the physics program. In any case, 
we anticipate that storage-ring designs are feasible that would allow 30-40% of the muon 
decays to provide useful neutrinos. 

Both Study-I and -II are site specific in that each has a few site-dependent aspects; 
otherwise, they are generic. In particular, Study-I assumed a new Fermilab booster to 
achieve its beam intensities and an underground storage ring. Study-II assumed BNL site- 
specific proton driver specifications corresponding to an upgrade of the 24-GeV AGS complex 
and a BNL-specific layout of the storage ring, which is housed in an above-ground berm to 
avoid penetrating the local water table. The primary substantive difference between the two 
studies is that Study-II aimed at a lower muon energy (20 GeV), but higher intensity (for 
physics reach) than Study-I. Taking the two Feasibility Studies together, we conclude that a 
high-performance Neutrino Factory could easily be sited at either BNL or Fermilab. Figure 
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^ shows a comparison of the performance of the Neutrino Factory designs in Study- 1 and 
Study-II j27| with the physics requirements. 
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FIG. 1: (Color) Muon decays in a straight section per 10 7 s vs. muon energy, with fluxes required 
for different physics searches assuming a 50 kT detector. Simulated performance of the two studies 
is indicated. 



To put the above performance figures in context, it is important to note that a /i + storage 
ring with an average neutrino energy of 15 GeV and 2 x 10 20 useful muon decays would yield 
(in the absence of oscillations) ~30,000 charged-current events in the v e channel per kiloton- 
year in a detector located 732 km away. In comparison, a 1.6 MW superbeam j3] from 
the Fermilab Main Injector with an average neutrino energy of 15 GeV would yield only 
~13,000 charged-current events per kiloton-year. In addition to having lower intensity 
than a Neutrino Factory beam, a superbeam would have significant v e contamination, which 
will be the major background in v ^ — > v e appearance searches. That is, it will be much easier 
to detect the oscillation v e — > from a muon storage ring neutrino beam than to detect 
the oscillation — > v e from a conventional neutrino beam, because the electron final state 
from the conventional beam has significant background contribution from 7r°'s produced in 



the events. 



C. Neutrino Factory Description 

The muons we use result from decays of pions produced when an intense proton beam 
bombards a high-power production target. The target and downstream transport channel 
are surrounded by superconducting solenoids to contain the pions and muons, which are pro- 
duced with a larger spread of transverse and longitudinal momenta than can be conveniently 
transported through an acceleration system. To prepare a beam suitable for subsequent ac- 
celeration, we first perform a "phase rotation," during which the initial large energy spread 
and small time spread are interchanged using induction linacs. Next, to reduce the trans- 
verse momentum spread, the resulting long bunch, with an average momentum of about 250 
MeV/c, is bunched into a 201.25-MHz bunch train and sent through an ionization cooling 
channel consisting of LH2 energy absorbers interspersed with rf cavities to replenish the 
energy lost in the absorbers. The resulting beam is then accelerated to its final energy using 
a superconducting linac to make the beam relativistic, followed by one or more recirculating 
linear accelerators (RLAs). Finally, the muons are stored in a racetrack-shaped ring with 
one long straight section aimed at a detector located at a distance of roughly 3000 km. A 
schematic layout is shown in Fig. |21 



D. Detector 

Specifications for the long-baseline Neutrino Factory detector are rather typical for an 
accelerator-based neutrino experiment. However, because of the need to maintain a high 
neutrino rate at these long distances (~3000 km), the detectors considered here are 3-10 
times more massive than those in current neutrino experiments. 

Several detector options could be considered for the far detector: 

• A 50 kton steel-scintillator-proportional-drift-tube (PDT) detector 

• A large water-Cherenkov detector, similar to SuperKamiokande but with either a 
magnetized water volume or toroids separating smaller water tanks j^]. 



A massive liquid-argon magnetized detector 
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FIG. 2: (Color) Schematic of the Neutrino Factory Study-II version. 

For the near detector, a compact liquid-argon TPC (similar to the ICARUS detector j^jj]) 
could be used. An experiment with a relatively thin Pb target (1 L ra( i), followed by a 
standard fixed-target spectrometer could also be considered. 

E. Staging Scenario 

If desired by the particle physics community, a fast-track plan leading directly to a Neu- 
trino Factory could be executed. On the other hand, the Neutrino Factory offers the distinct 
advantage that it can be built in stages. This could satisfy both programmatic and cost 
constraints by allowing an ongoing physics program while reducing the annual construction 
funding needs. Depending on the results of our technical studies and the results of ongoing 
searches for the Higgs boson, it is hoped that the Neutrino Factory is really the penultimate 
stage, to be followed later by a Muon Collider (e.g., a Higgs Factory). Such a collider offers 
the potential of bringing the energy frontier in particle physics within reach of a moderate- 
sized machine. Possible stages for the evolution of a muon beam facility are described in 
Section rill .11 
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F. R&D Program 



Successful construction of a muon storage ring to provide a copious source of neutrinos 
requires development of many novel approaches; construction of a high-luminosity Muon 
Collider requires even more. It was clear from the outset that the breadth of R&D issues 
to be dealt with would be beyond the resources available at any single national laboratory 
or university For this reason, in 1995, interested members of the high-energy physics and 
accelerator physics communities formed the MC to coordinate the required R&D efforts 
nationally. The task of the MC is to define and carry out R&D needed to assess the technical 
feasibility of constructing initially a muon storage ring that will provide intense neutrino 
beams aimed at detectors located many thousands of kilometers from the accelerator site, 
and ultimately a fi + n~ collider that will carry out fundamental experiments at the energy 
frontier in high-energy physics. 

The MC also serves to coordinate muon-related R&D activities of the NSF-sponsored 
University Consortium (UC) and the state-sponsored Illinois Consortium for Accelerator 
Research (ICAR), and is the focal point for defining the needs of muon-related R&D to 
the managements of the sponsoring national laboratories and to the funding agencies (both 
DOE and NSF). As already noted, though the MC was formed initially to carry out R&D 
that might lead eventually to the construction of a Muon Collider, more recently its focus 
has shifted mainly, but not exclusively, to a Neutrino Factory. 

The MC maintains close contact with parallel R&D efforts under way in Europe (cen- 
tered at CERN) and in Japan (centered at KEK). Through its international members, the 
MC also fosters coordination of the international muon-beam R&D effort. Two major ini- 
tiatives, a Targetry Experiment (E951) in operation at BNL and a Muon Cooling R&D 
program (MUCOOL), have been launched by the MC. In addition, the Collaboration, work- 
ing in conjunction with the UC and ICAR in some areas, coordinates substantial efforts in 
accelerator physics and component R&D to define and assess parameters for feasible designs 
of muon-beam facilities. 
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G. Outline of Report 



In what follows, we give the motivation and a scenario for a staged approach to con- 
structing a Neutrino Factory and eventually a Muon Collider. Section |H] discusses the 
physics opportunities, starting from conventional "superbeams" and going to cold muon 
beams, then a Neutrino Factory with its near and far detectors, and finally a Muon Collider. 
In Section IIII1 we describe the components of a Neutrino Factory, based on the Study-II 
design, and indicate a scientifically productive staged path for reaching it. Section ITVl cov- 
ers our present concept of an entry-level Higgs Factory Muon Collider. In support of the 
construction of a Neutrino Factory, an R&D program is already under way to address vari- 
ous technical issues. A description of the status and plans for this program is presented in 
Section Section IVT1 describes current thinking about a cooling demonstration experiment 
that would be carried out as an international effort. Finally, in Section fVIII we provide a 
brief summary of our work. 



II. PHYSICS MOTIVATION 



In this Section we cover the physics potential of the Neutrino Factory accelerator complex, 
which includes superbeams of conventional neutrinos that are possible using the proton 
driver needed for the factory, and intense beams of cold muons that become available once 
the muon cooling and collection systems for the factory are in place. Once the cold muons 
are accelerated and stored in the muon storage ring, we realize the full potential of the 
factory in both neutrino oscillation and non-oscillation physics. Cooling muons will be a 
learning experience. We hope that the knowledge gained in constructing a Neutrino Factory 
can be used to cool muons sufficiently to produce the first muon collider operating as a 
Higgs factory. We examine the physics capabilities of such a collider, which if realized, will 
invariably lead to higher energy muon colliders with exciting physics opportunities. 



A. Neutrino Oscillation Physics 

Here we discuss (33 1 the current evidence for neutrino oscillations, and hence neutrino 
masses and lepton mixing, from solar and atmospheric data. A review is given of some theo- 
retical background including models for neutrino masses and relevant formulas for neutrino 
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oscillation transitions. We next mention the near-term and mid-term experiments in this 
area and comment on what they hope to measure. We then discuss the physics potential of 
a muon storage ring as a Neutrino Factory in the long term. 



1. Evidence for Neutrino Oscillations 

In a modern theoretical context, one generally expects nonzero neutrino masses and asso- 
ciated lepton mixing. Experimentally, there has been accumulating evidence for such masses 
and mixing. All solar neutrino experiments (Homestake, Kamiokande, SuperKamiokande, 
SAGE, GALLEX and SNO) show a significant deficit in the neutrino fluxes coming from the 



Sun [34j. This deficit can be explained by oscillations of the i/ e 's into other weak eigenstate(s), 
with Ami, of the order 10~ 5 eV 2 for solutions involving the Mikheyev-Smirnov-Wolfenstein 

nn 

(MSW) resonant matter oscillations |3j| |3J| or of the order of 10 10 eV 2 for vacuum oscilla- 
tions [3J|. Accounting for the data with vacuum oscillations (VO) requires almost maximal 
mixing. The MSW solutions include one for small mixing angle (SMA) and one for large 
mixing angle (LMA). Another piece of evidence for neutrino oscillations is the atmospheric 



neutrino anomaly, observed by Kamiokande j^J, 1MB 41 
highest statistics, and by Soudan Q and MACRO 



SuperKamiokande 



4J with the 



44j . These data can be fit by the 



inference of — > v x oscillations with Am 2 tm ~ 3 x 10~ 3 eV 2 |42| and maximal mixing 
sin 2 26' a t m = 1. The identification v x = v T is preferred over v x = ^sterile, and the identifica- 
tion v x = v e is excluded by both the Superkamiokande data and the Chooz experiment 4^j] . 
In addition, the LSND experiment has reported — > v e and — > u e oscillations with 
Am 2 SND ~ 0.1-1 eV 2 and a range of possible mixing angles. This result is not confirmed, but 
also not completely ruled out, by a similar experiment, KARMEN J^tJ. The miniBOONE 
experiment at Fermilab is designed to resolve this issue, as discussed below. If one were to 
try to fit all of these experiments, then, since they involve three quite different values of 
Amfj = m(ui) 2 — m(uj) 2 , which could not satisfy the identity for three neutrino species, 

Aml 2 + Am^ + Am 2 3 = , (1) 

it would follow that one would have to introduce at least one further neutrino. Since it is 
known from the measurement of the Z width that there are only three leptonic weak doublets 
with associated light neutrinos, it follows that such further neutrino weak eigenstate(s) would 
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have to be electroweak singlet(s) ("sterile" neutrinos). Because the LSND experiment has 
not been confirmed by the KARMEN experiment, we choose here to use only the (confirmed) 
solar and atmospheric neutrino data in our analysis, and hence to work in the context of 
three active neutrino weak eigenstates. 



2. Neutrino Oscillation Formalism 



In this theoretical context, consistent with solar and atmospheric data, there are three 
electroweak-doublet neutrinos and the neutrino mixing matrix is described by 



C12C13 C13S12 si 3 e lS ' 



U 



— C23S12 — Si3S23Ci2e J<5 C 12C23 — s 12 s 13 s 23 el<5 C13S23 
S12-S23 — •5i3C 12 C23e* <5 — S23C 12 — S^C^S^e** 5 C13C23 



K> , (2) 



/ 



where Qj = cos#y, = sin#y, and K' = diag(l, e 1 ^ 1 , e 1 ^ 2 ). The phases 0i and 02 do 
not affect neutrino oscillations. Thus, in this framework, the neutrino mixing relevant for 
neutrino oscillations depends on the four angles 6 , 12 , 6 , 13 , # 23 , and 5, and on two independent 
differences of squared masses, Am 2 tm , which is Am| 2 = m{v^) 2 — m(i/ 2 ) 2 in the favored 
fit, and Am^ ol , which may be taken to be Am 21 = m(z/ 2 ) 2 — m(ui) 2 . Note that these Am 2 
quantities involve both magnitude and sign; although in a two-species neutrino oscillation in 
vacuum the sign does not enter, in the three-species-oscillation, which includes both matter 
effects and CP violation, the signs of the Am 2 quantities enter and can, in principle, be 
measured. For our later discussion it will be useful to record the formulas for the various 
neutrino-oscillation transitions. In the absence of any matter effect, the probability that a 
(relativistic) weak neutrino eigenstate u a becomes v b after propagating a distance L is 

-Am?L- 



P{v a -> u b ) = 5 ab -iYl Re (Kab,ij) sin^^^J + 

i>j=l 

/ Am? - L\ / Am?- L \ 
+ 4 ^ M#«w) sin(^) cos(^) (3) 



i>j=l 



where 



and 



K ab , l3 = U ai U* K U* aj U bj (4) 



Am 2 - = m(vi) 2 — m{vjf . (5) 
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Recall that in vacuum, CPT invariance implies P{pb v a ) = P[y a — > Vb) and hence, for 
b = a, P{p a v a ) = P(v a — > u a ). For the CP-transformed reaction v a — > Vb and the 
T-reversed reaction — > u a , the transition probabilities are given by the right-hand side of 
(J2J) with the sign of the imaginary term reversed. (Below we shall assume CPT invariance, 
so that CP violation is equivalent to T violation.) In most cases there is only one mass scale 
relevant for long-baseline neutrino oscillations, Am 2 tm ~ few x 10~ 3 eV 2 , and one possible 
neutrino mass spectrum is the hierarchical one 

Am 2 21 = Am 2 sol < Am^ « Am 2 2 = Am^ tm . (6) 

In this case, CP (T) violation effects may be negligibly small, so that in vacuum 

P(Va -> h) = P(Va -> V b ) (7) 

and 

P(y b -> Va) = P{Va -> Vb) ■ (8) 

In the absence of T violation, the second equality (jSJ) would still hold in uniform matter, 
but even in the absence of CP violation, the first equality (JJJ) would not hold. With the 
hierarchy the expressions for the specific oscillation transitions are 

P(u t 



PiVe 



In neutrino oscillation searches using reactor antineutrinos, i.e, tests of v e — > u e , the two- 
species mixing hypothesis used to fit the data is 

P[v e ->V e ) = P("e -> V X ) 

X 

= 1 - sin 2 (2£ rcactor ) sin 2 ( Am ~ L ) , (12) 
1(3 



V-r) = 4|t/33| 2 |^3| 2 sin 2 (^|^) 

= sin 2 (2^ 23 ) cos 4 (M sin 2 (^f^) , (9) 

v,) = 4|^ 3 | 2 |^3| 2 sin 2 (^|^) 

= sin 2 (2^ 13 ) sin 2 (£ 23 ) sin 2 (^|^) , (10) 

*r) = 4|f/33| 2 |f/l 3 | 2 sin 2 (^|^) 

= sin 2 (2^ 3 )cos 2 (^3)sin 2 (^|^) . (11) 



where Am 2 cactor is the squared mass difference relevant for v e — > u x . In particular, in the 
upper range of values oi 
disappearance, one has 



upper range of values of Am atm , since the transitions v e — > z/„ and v e — > v T contribute to u, 



P(v e - v e ) = 1 - sin 2 (2£ 13 ) sin 2 (^f^ ) , (13) 
i.e., Reactor = #13, and, for the value |Am| 2 | = 3 x 10_^ 3 eV 2 from SuperK, the CHOOZ 



experiment on v e disappearance yields the upper limit 

sin 2 (2# 13 ) <0.1, (14) 



which is also consistent with conclusions from the SuperK data analysis 42j. Further, the 
quantity "sin 2 (2# atm )" often used to fit the data on atmospheric neutrinos with a simplified 
two-species mixing hypothesis, is, in the three-generation case, 

sin 2 (2# atm ) = sin 2 (2# 23 ) cos 4 (# 13 ) . (15) 

The SuperK experiment finds that the best fit to their data is — > v T oscillations with 
maximal mixing, and hence sin 2 (2^23) = 1 and |#i 3 | <C 1. The various solutions of the solar 
neutrino problem involve quite different values of Am^ and sin 2 (26 l i2): (i) large mixing angle 
solution, LMA: Am| x ~ few x 10~ 5 eV 2 and sin 2 (2#i2) — 0.8; (ii) small mixing angle solution, 
SMA: Am^ ~ 10" 5 eV 2 and sin 2 (2^ 2 ) ~ 10" 2 , (iii) LOW: Am 2 21 ~ 10" 7 eV 2 , sin 2 (2^ 2 ) ~ 1, 
and (iv) "just-so": Am 21 ~ 10~ 10 eV 2 , sin 2 (2# 12 ) ~ 1. The SuperK experiment favors the 
LMA solutions [34]]; for other global fits, see, e.g., Ref. (34)]. We have reviewed the three 
neutrino oscillation phenomenology that is consistent with solar and atmospheric neutrino 
oscillations. In what follows, we will examine the neutrino experiments planned for the 
immediate future that will address some of the relevant physics. We will then review the 
physics potential of the Neutrino Factory. 



3. Relevant Near- and Mid-Term Experiments 

There are currently intense efforts to confirm and extend the evidence for neutrino os- 
cillations in all of the various sectors — solar, atmospheric, and accelerator. Some of these 
experiments are running; in addition to SuperKamiokande and Soudan-2, these include the 
Sudbury Neutrino Observatory, SNO, and the K2K long baseline experiment between KEK 
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and Kamioka. Others are in development and testing phases, such as miniBOONE, MINOS, 
the CERN-Gran Sasso program, KamLAND, Borexino, and MONOLITH [48]. Among the 
long baseline neutrino oscillation experiments, the approximate distances are L ~ 250 km for 
K2K, 730 km for both MINOS (from Fermilab to Soudan) and the proposed CERN-Gran 
Sasso experiments. K2K is a disappearence experiment with a conventional neutrino 
beam having a mean energy of about 1.4 GeV, going from KEK 250 km to the SuperK 
detector. It has a near detector for beam calibration. It has obtained results consistent 
with the SuperK experiment, and has reported that its data disagree by 2a with the no- 
oscillation hypothesis [3]. MINOS is another conventional neutrino beam experiment that 
takes a beam from Fermilab 730 km to a detector in the Soudan mine in Minnesota. It 
again uses a near detector for beam flux measurements and has opted for a low-energy con- 
figuration, with the flux peaking at about 3 GeV. This experiment is scheduled to start 
taking data in 2005 and, after some years of running, to obtain higher statistics than the 
K2K experiment and to achieve a sensitivity down to the level [Amiy ~ 10~ 3 eV 2 . The 
CERN-Gran Sasso program will also come on in 2005. It will use a higher-energy neutrino 
beam, E v ~ 17 GeV, from CERN to the Gran Sasso deep underground laboratory in Italy. 
This program will emphasize detection of the r's produced by the z/ r 's that result from the 
inferred neutrino oscillation transition — > v T . The OPERA experiment will do this using 
emulsions j^j], while the ICARUS proposal uses a liquid argon chamber j^. For the joint 
capabilities of MINOS, ICARUS and OPERA experiments see Ref. Plans for the Japan 
Hadron Facility (JHF), also called the High Intensity Proton Accelerator (HIPA), include 
the use of a 0.77 MW proton driver to produce a high-intensity conventional neutrino beam 
with a path length of 300 km to the SuperK detector (5^. Moreover, at Fermilab, the mini- 
BOONE experiment is scheduled to start data taking in the near future and to confirm or 
refute the LSND claim after a few years of running. There are several neutrino experiments 
relevant to the solar neutrino anomaly. The SNO experiment is currently running and has 
recently reported their first results that confirm solar neutrino oscillations . These involve 
measurement of the solar neutrino flux and energy distribution using the charged current 
reaction on heavy water, u e + d — > e + p + p. They are expected to report on the neutral 
current reaction v e + d — > v e + n + p shortly. The neutral current rate is unchanged in 
the presence of oscillations that involve standard model neutrinos, since the neutral current 
channel is equally sensitive to all the three neutrino species. If however, sterile neutrinos are 
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involved, one expects to see a depletion in the neutral current channel also. However, the un- 
certain normalization of the 8 B flux makes it difficult to constrain a possible sterile neutrino 
component in the orations Q. The KamLAND experiment Q in Japan 8 ta rt ed taking 
data in January 2002. This is a reactor antineutrino experiment using baselines of 100-250 
km. It will search for v e disappearance and is sensitive to the solar neutrino oscillation scale. 
KamLAND can provide precise measurements of the LMA solar parameters and recently 
the first results from KamLAND have confirmed the LMA solution js^j. A global analysis 
of the KamLAND and solar neutrino data has further restricted the solar b~m? range and 
the best ht vaine on.ently is 7 X !„- eV* HQ. On a similar tnne s ca,e, the Borexano 
experiment in Gran Sasso is scheduled to turn on and measure the 7 Be neutrinos from the 
sun. These experiments should help us determine which of the various solutions to the solar 
neutrino problem is preferred, and hence the corresponding values of Am^ and sin 2 (2(9i2). 
This, then, is the program of relevant experiments during the period 2000-2010. By the end 
of this period, we may expect that much will be learned about neutrino masses and mixing. 
However, there will remain several quantities that will not be well measured and which can 
be measured by a Neutrino Factory. 



4. Oscillation Experiments at a Neutrino Factory 

Although a Neutrino Factory based on a muon storage ring will turn on several years 
after this near-term period in which K2K, MINOS, and the CERN-Gran Sasso experiments 
will run, it has a valuable role to play, given the very high-intensity neutrino beams of fixed 
flavor-pure content, including, uniquely, v e and v e beams in addition to and beams. 
A conventional positive charge selected neutrino beam is primarily with some admixture 
of iVs and other flavors from K decays (0(1%) of the total charged current rate) and the 
fluxes of these neutrinos can only be fully understood after measuring the charged particle 
spectra from the target with high accuracy. In contrast, the potential of the neutrino beams 
from a muon storage ring is that the neutrino beams would be of extremely high purity: 
yU~ beams would yield 50% and 50% v e , and /i + beams, the charge conjugate neutrino 
beams. Furthermore, these could be produced with high intensities and low divergence that 
make it possible to go to longer baselines. 

In what follows, we shall take the design values from Study-II of 10 20 fi decays per 
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TABLE I: Neutrino-oscillation modes that can be studied with conventional neutrino beams or 
with beams from a Neutrino Factory, with ratings as to degree of difficulty in each case; * = well 
or easily measured, 1/ = measured poorly or with difficulty, — = not measured. 
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"Snowmass year" (10 7 sec) as being typical. The types of neutrino oscillations that can be 
searched for with the Neutrino Factory based on the muon storage ring are listed in Table H] 
for the case of \x~ which decays to z/ M e~z/ e : It is clear from the processes listed that since the 
beam contains both neutrinos and antineutrinos, the only way to determine the flavor of the 
parent neutrino is to determine the identity of the final state charged lepton and measure 
its charge. A capability unique to the Neutrino Factory will be the measurement of the 
oscillation v e — > u^, giving a wrong-sign Of greater difficulty would be the measurement 
of the transition v e — > v T , giving a r + which will decay part of the time to fi + . These 
physics goals mean that a detector must have excellent capability to identify muons and 
measure their charges. Especially in a steei-scintillator detector, the oscillation — > v e 
would be difficult to observe, since it would be difficult to distinguish an electron shower 
from a hadron shower. From the above formulas for oscillations, one can see that, given the 
knowledge of |Am| 2 | and sin 2 (2# 2 3) that will be available by the time a Neutrino Factory 
is built, the measurement of the v e — > transition yields the value of $13. To get a rough 
idea of how the sensitivity of an oscillation experiment would scale with energy and baseline 
length, recall that the event rate in the absence of oscillations is simply the neutrino flux 
times the cross section. First of all, neutrino cross sections in the region above about 10 GeV 
(and slightly higher for r production) grow linearly with the neutrino energy. Secondly, the 
beam divergence is a function of the initial muon storage ring energy; this divergence yields 

20 



a flux, as a function of 64, the angle of deviation from the forward direction, that goes like 
~ E 2 . Combining this with the linear E dependence of the neutrino cross section and 
the overall 1/L 2 dependence of the flux far from the production region, one finds that the 
event rate goes like 

We base our discussion on the event rates given in the Fermilab Neutrino Factory study (27 1 . 
For a stored muon energy of 20 GeV, and a distance of L = 2900 to the WIPP Carlsbad site 
in New Mexico, these event rates amount to several thousand events per kton of detector per 
year, i.e, they are satisfactory for the physics program. This is also true for the other path 
lengths under consideration, namely L = 2500 km from BNL to Homestake and L = 1700 km 



ow a single pathlength L, but a bowtie 
6(J). We anticipate that at a time when 



to Soudan. A usual racetrack design would only ah 
design could allow two different path lengths (e.g., 
the Neutrino Factory turns on, |Am 2 2 | and sin 2 (2# 23 ) would be known at perhaps the 10% 
level (while recognizing that future projections such as this are obviously uncertain). The 
Neutrino Factory will significantly improve precision in these parameters, as can be seen from 
Fig. 01 which shows the error ellipses possible for a 30 GeV muon storage ring. In addition, 
the Neutrino Factory can contribute to the measurement of: (i) #13 , as discussed above; (ii) 
measurement of the sign of Am^ using matter effects; and (iii) possibly a measurement of 
CP violation in the leptonic sector, if sin 2 (2# 13 ), sin 2 (2# 21 ), and Am| x are sufficiently large. 
To measure the sign of Am 2 2 , one uses the fact that matter effects reverse sign when one 
switches from neutrinos to antineutrinos, and carries out this switch in the charges of the 
stored n ± . We elaborate on this next. 



5. Matter Effects 

With the advent of the muon storage ring, the distances at which one can place detectors 
are large enough so that for the first time matter effects can be exploited in accelerator-based 
oscillation experiments. Simply put, matter effects are the matter- induced oscillations that 
neutrinos undergo along their flight path through the Earth from the source to the detector. 
Given the typical density of the earth, matter effects are important for the neutrino energy 
range E ~ O(10) GeV and Am§ 2 ~ 10~ 3 eV 2 , values relevant for the long baseline exper- 
iments. Matter effects in neutrino propagation were first pointed out by Wolfenstein [3] 
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FIG. 3: Fit to muon neutrino survival distribution for = 30 GeV and L = 2800 km for 10 pairs 
of sin 2 2#, <5m 2 values. For each fit, the lcr, 2a and 3a contours are shown. The generated points 
are indicated by the dark rectangles and the fitted values by stars. The SuperK 68%, 90%, and 
99% confidence levels are superimposed. Each point is labelled by the predicted number of signal 
events for that point. 

and Barger, Pakvasa, Phillips and Whisnant j^]. (See the papers (fill jzfl for details of 
the matter effects and their relevance to neutrino factories.) In brief, assuming a normal 
hierarchy, the transition probabilities for propagation through matter of constant density 
are HQ 



P(is„ 



2„2 



x f + 2xyfg(cos 5 cos A + sin 5 sin A) + y g 



2„2 



cot 623X f — 2xyfg(cos5 cos A + sin <5sin A) + tan Q^V 9 
sin 2 26*23 sin 2 A 



(17) 
(18) 
(19) 



+ a sin 26*23 sin 2 A ( sin 6*13 sin 2#i2 cos 26*23 sin A — A cos 2 6u sin 2^23 | 
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where 
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\Sm 2 31 \L/4:E u = 1.27|<5m2 1 /eV 2 |(L/km)/(^/GeV) , 
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(21) 


CM 


IAtT7^ 1 F\TY^ 1 

l°'"'2l/ w '' t 3ll ' 




X EE 


sin #23 sin 2#i 3 , 


(23) 
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a cos #23 sin 2#i 2 , 


(24) 


f = 


sin((l T A)A)/(l T i), 


(25) 


9 = 


sin(iA)/i. 


(26) 



The amplitude A for v e e forward scattering in matter is given by 

A = 2V2G F N e E u = 1.52 x 1(T 4 eV 2 Y e p( g/cm 3 )£( GeV) . (27) 

Here Y e is the electron fraction and p(x) is the matter density. For neutrino trajectories 
that pass through the earth's crust, the average density is typically of order 3 gm/cm 3 and 
Y e ~ 0.5. For neutrinos with Sm^ > or anti-neutrinos with Sm^ < 0, A = 1 corresponds 
to a matter resonance. Thus, for a Neutrino Factory operating with positive stored muons 
(producing a u e beam) one expects an enhanced production of opposite sign (p~) charged- 
current events as a result of the oscillation v e —> if 5m| 2 is positive and vice versa for 
stored negative beams. Figure 0] jO] shows the wrong-sign muon appearance spectra as 
function of Sm 2 2 for both p + and p~ beams for both signs of 8m\ 2 at a baseline of 2800 km. 
The resonance enhancement in wrong sign muon production is clearly seen in Fig. Efb) and 
(c). By comparing these (using first a stored p + beam and then a stored p~ beam) one 
can thus determine the sign of Am| 2 as well as the value of sin 2 (2#i3). Figure El [3] shows 
the difference in negative log-likelihood between a correct and wrong-sign mass hypothesis 
expressed as a number of equivalent Gaussian standard deviations versus baseline length for 
muon storage ring energies of 20, 30, 40 and 50 GeV. The values of the oscillation parameters 
are for the LMA scenario with sin 2 29 13 = 0.04. Figure |3fa) is for 10 20 decays for each sign 
of stored energy and a 50 kiloton detector and positive 5m\ 2 , (b) is for negative 8m 2 2 for 
various values of stored muon energy. Figures El (c) and (d) show the corresponding curves 
for 10 19 decays and a 50 kiloton detector. An entry-level machine would permit one to 
perform a 5a differentiation of the sign of 8m\ 2 at a baseline length of ~2800 km. For the 
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FIG. 4: (Color)The wrong sign muon appearance rates for a 20 GeV muon storage ring at a baseline 
of 2800 km with 10 20 decays and a 50 kiloton detector for (a) /U + stored and negative 6m^ 2 > (b) [i~ 
stored and negative b~m\ 2 , (c) ^ + stored and positive 5m| 2 , (d) \i~ stored and positive bm\ 2 . The 
values of |<5m| 2 | range from 0.0005 to 0.0050 eV 2 in steps of 0.0005 eV 2 . Matter enhancements are 
evident in (b) and (c). 

Study II design, in accordance with the previous Fermilab study [^tJ, one estimates that it 
is possible to determine the sign of 8m\ 2 even if sin 2 (26'i3) is as small as ~ 1CT 3 . 

6. CP Violation 

CP violation is measured by the (rephasing-invariant) product 
J = Im(U ai U^U* aj U bj ) 

= - sin(2# 12 ) sin(2# 13 ) cos(^ 3 ) sin(2# 23 ) sin 5 . (28) 
8 

Leptonic CP violation also requires that each of the leptons in each charge sector be non- 
degenerate with any other leptons in this sector; this is, course, true of the charged lepton 
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FIG. 5: The statistical significance (number of standard deviations) with which the sign of <5mJ 2 
can be determined versus baseline length for various muon storage ring energies. The results are 
shown for a 50 kiloton detector, and (a) 10 20 // + and [i~ decays and positive values of 5m| 2 ; (b) 10 20 
fi + and n~ decays and negative values of 5m| 2 ; (c) 10 19 /U + and ff~ decays and positive values of 
(5 m 32) (d) 10 19 /i + and fx~ decays and negative values of <5?n 2 2 . 

sector and, for the neutrinos, this requires Arnf,- ^ for each such pair ij. In the quark 
sector, J is known to be small: Jckm ~ O(10~ 5 ). A promising asymmetry to measure is 
P[y e — > v^) — P{y e — v^). As an illustration, in the absence of matter effects, 

P[y e -=> Vy) - P[p e -=> = 4J(sin20 32 + sin20 2 i + sin20i 3 ) 

= -16Jsin0 32 sin0i 3 sin0 2 i , (29) 
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where 



Am?L 



(30) 



AE 



In order for the CP violation in Eq. fj29jl to be large enough to measure, it is necessary 
that #12, 6*i3, and Am s 2 | = Am^ not be too small. From atmospheric neutrino data, we 
have 023 — vr/4 and #13 <C 1. If LMA describes solar neutrino data, then sin 2 (2^12) — 0.8, 
so J ~ 0.1 sin(26* 13 ) sin 5. For example, if sin 2 (2# 13 ) = 0.04, then J could be ^> Jckm- 
Furthermore, for parts of the LMA phase space where Am 2 ol ~ 4x 10~ 5 eV 2 the CP violating 
effects might be observable. In the absence of matter, one would measure the asymmetry 



However, in order to optimize this ratio, because of the smallness of Am 2 x even for the LMA, 
one must go to large pathlengths L, and here matter effects are important. These make 
leptonic CP violation challenging to measure, because, even in the absence of any intrinsic 
CP violation, these matter effects render the rates for v e — > and v e — > unequal since 
the matter interaction is opposite in sign for v and v. One must therefore subtract out 
the matter effects in order to try to isolate the intrinsic CP violation. Alternatively, one 
might think of comparing v e — > v ^ with the time-reversed reaction — > v e . Although this 
would be equivalent if CPT is valid, as we assume, and although uniform matter effects are 
the same here, the detector response is quite different and, in particular, it is quite difficult 
to identify e^. Results from SNO and KamLAND testing the LMA |5J| will help further 
planning. The Neutrino Factory provides an ideal set of controls to measure CP violation 
effects since we can fill the storage ring with either /i + or /i~ particles and measure the 
ratio of the number of events v e — > v^jve *V Figure |H1 shows this ratio for a Neutrino 
Factory with 10 21 decays and a 50 kiloton detector as a function of the baseline length. 
The ratio depends on the sign of <5m 2 2 . The shaded band around either curve shows the 
variation of this ratio as a function of the CP- violating phase 5. The number of decays 
needed to produce the error bars shown is directly proportional to sin 2 #13, which for the 
present example is set to 0.004. Depending on the magnitude of J, one may be driven to 
build a Neutrino Factory just to understand CP violation in the lepton sector, which could 
have a significant role in explaining the baryon asymmetry of the Universe |78j . 
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FIG. 6: (Color)Predicted ratios of wrong-sign muon event rates when positive and negative muons 
are stored in a 20 GeV Neutrino Factory, shown as a function of baseline. A muon measurement 
threshold of 4 GeV is assumed. The lower and upper bands correspond, respectively, to negatve 
and positive Sm^- The widths of the bands show how the predictions vary as the CP violating 
phase 5 is varied from —tt/2 to tt/2, with the thick lines showing the predictions for 5 = 0. The 
statistical error bars correspond to a high-performance Neutrino Factory yielding a data sample of 
10 21 decays with a 50 kiloton detector. The curves are based on calculations presented in [721 ] . 

B. Physics Potential of Superbeams 

It is possible to extend the reach of the current conventional neutrino experiments by 
enhancing the capabilities of the proton sources that drive them. These enhanced neutrino 
beams have been termed "superbeams" and form an intermediate step on the way to a 



Neutrino Factory. Their capabilities have been explored in recent papers (28 



79 



80JI. These 



articles consider the capabilities of enhanced proton drivers at (i) the proposed 0.77 MW 



50 GeV proton synchrotron at the Japan Hadron Facility (JHF) 



531 ] , (ii) a 4 MW upgraded 
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version of the JHF, (iii) a new ~ 1 MW 16 GeV proton driver |8l( that would replace the 
existing 8 GeV Booster at Fermilab, or (iv) a fourfold intensity upgrade of the 120 GeV Fer- 
milab Main Injector (MI) beam (to 1.6 MW) that would become possible once the upgraded 
(16 GeV) Booster was operational. Note that the 4 MW 50 GeV JHF and the 16 GeV 
upgraded Fermilab Booster are both suitable proton drivers for a neutrino factory. The 
conclusions of both reports are that superbeams will extend the reaches in the oscillation 
parameters of the current neutrino experiments but "the sensitivity at a Neutrino Factory to 
CP violation and the neutrino mass hierarchy extends to values of the amplitude parameter 
sin 2 2^13 that are one to two orders of magnitude lower than at a superbeam" 79. 8o|. 



To illustrate these points, we choose one of the most favorable superbeam scenarios 
studied: a 1.6 MW NuMI-like high energy beam with L = 2900 km, detector param- 
eters corresponding to the liquid argon scenario in [79! . IsO^ . and oscillation parameters 
l^ m 32l = 3-5 x 10 -3 eV 2 and Sm^ = 1 x 10~ 4 eV 2 . The calculated three-sigma error 
ellipses in the (N(e + ), N(e~)) plane are shown in Fig. for both signs of 5m 2 2 , with the 
curves corresponding to various CP phases 5 (as labeled). The magnitude of the — > v e 
oscillation amplitude parameter sin 2 26\z varies along each curve, as indicated. The two 
groups of curves, which correspond to the two signs of 5m 2 2 , are separated by more than 
3cr provided sin 2 2^13 > 0.01. Hence the mass heirarchy can be determined provided the 
v ii v e oscillation amplitude is not more than an order of magnitude below the currently 
excluded region. Unfortunately, within each group of curves, the CP-conserving predictions 
are separated from the maximal CP-violating predictions by at most 3a. Hence, it will be 
difficult to conclusively establish CP violation in this scenario. Note for comparison that a 
very long baseline experiment at a neutrino factory would be able to observe v e —* os- 
cillations and determine the sign of Sm^ 2 for values of sin 2 2^13 as small as (9(0.0001). This 
is illustrated in Fig. |HJ A Neutrino Factory thus outperforms a conventional superbeam in 
its ability to determine the sign of 8m\ 2 - Comparing Fig. and Fig. |H1 one sees that the 
value of sin 2 2^13, which has yet to be measured, will determine the parameters of the first 
Neutrino Factory. 

Finally, we compare the superbeam — * v e reach with the corresponding Neutrino Fac- 
tory v e — > reach in Fig. El which shows the 3<r sensitivity contours in the (5m 2 !, sin 2 2# 13 ) 
plane. The superbeam sin 2 26 , 13 reach of a few x 10 -3 is almost independent of the sub-leading 
scale Sm^i- However, since the neutrino factory probes oscillation amplitudes O(10 -4 ) the 
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FIG. 7: Three-sigma error ellipses in the (N(e + ), N(e~)) plane, shown for — > u e and P„ — > v e 
oscillations in a NuMI-like high energy neutrino beam driven by a 1.6 MW proton driver. The 
calculation assumes a liquid argon detector with the parameters listed in [2^ . a baseline of 2900 km, 
and 3 years of running with neutrinos, 6 years running with antineutrinos. Curves are shown for 
different CP phases 5 (as labelled), and for both signs of <5m| 2 with |5m| 2 | = 0.0035 eV 2 , and the 
sub-leading scale Sm^i = 10 -4 eV 2 . Note that sin 2 2#i3 varies along the curves from 0.001 to 0.1, 
as indicated [^J. 

sub-leading effects cannot be ignored, and v e —>■ events would be observed at a Neutrino 
Factory over a significant range of Srn^i even if sin 2 2# 13 = 0. 



C. Non-oscillation physics at a Neutrino Factory 



The study of the utility of intense neutrino beams from a muon stora ge ri ng in determining 
the parameters governin g n on-oscillation physics was begun in 1997 14j . More complete 



studies can be found in 



and recently a European group has brought out an extensive 



study on this topic |82j. A Neutrino Factory can measure individual parton distributions 
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FIG. 8: Three-sigma error ellipses in the (iV(//+), iV(/x— )) plane, shown for a 20 GeV neutrino 
factory delivering 3.6 x 10 21 useful muon decays and 1.8 x 10 21 antimuon decays, with a 50 kt 
detector at L = 7300 km, drn^i = 10 -4 eV 2 , and 5 = 0. Curves are shown for both signs of <5m 2 2 ; 
sin 2 20x3 varies along the curves from 0.0001 to 0.01, as indicated [7^ . 

within the proton for all light quarks and anti-quarks. It could improve valence distributions 
by an order of magnitude in the kinematical range x > 0.1 in the unpolarized case. The 
individual components of the sea (u, d, s and s) , as well as the gluon, would be measured 
with relative accuracies in the range of 1-10%, for 0.1 < x < 0.6. A full exploitation of the 
Neutrino Factory potential for polarized measurements of the shapes of individual partonic 
densities requires an a priori knowledge of the polarized gluon density. The forthcoming set 
of polarized deep inelastic scattering experiments at CERN, DESY and RHIC may provide 
this information. The situation is also very bright for measurements of C-even distributions. 
Here, the first moments of singlet, triplet and octet axial charges can be measured with 
accuracies that are up to one order of magnitude better than the current uncertainties. In 
particular, the improvement in the determination of the singlet axial charge would allow a 
definitive confirmation or refutation of the anomaly scenario compared to the 'instanton' 
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FIG. 9: Summary of the 3d level sensitivities for the observation of Vn — > v e at various MW-scale 
superbeams (as indicated) with liquid argon "A" and water cerenkov "W" detector parameters, 
and the observation of v e — > in a 50 kt detector at 20, 30, 40, and 50 GeV neutrino factories 
delivering 2 x 10 20 muon decays in the beam-forming straight section. The limiting 3<7 contours 
are shown in the (Sm^i, sin 2 2613) plane. All curves correspond to 3 years of running. The grey 
shaded area is already excluded by current experiments. 

or 'skyrmion' scenarios, at least if the theoretical uncertainty originating from the small- x 
extrapolation can be kept under control. The measurement of the octet axial charge with a 
few percent uncertainty will allow a determination of the strange contribution to the proton 
spin better than 10%, and allow stringent tests of models of SU(3) violation when compared 
to the direct determination from hyperon decays. A measurement of as(Mz) and sin 2 9w 
will involve different systematics from current measurements and will therefore provide an 
important consistency check of current data, although the accuracy of these values is not 
expected to be improved. The weak mixing angle can be measured in both the hadronic 
and leptonic modes with a precision of approximately 2 x 10~ 4 , dominated by the statistics 



Excluded 



0.77 MW JHF , W , L = 295 km 



4 MW SJHF , A , L = 295 km 



1 .6 MW SNuMI ME , A , L = 2900 km 
NuFact L = 2900 km 




31 



and the luminosity measurement. This determination would be sensitive to different classes 
of new-physics contributions. Neutrino interactions are a very good source of clean, sign- 
tagged charm particles. A Neutrino Factory can measure charm production with raw event 
rates up to 100 million charm events per year with ~ 2 million double-tagged events. (Note 
that charm production becomes significant for storage ring energies above 20 GeV). Such 
large samples are suitable for precise extractions of branching ratios and decay constants, 
the study of spin-transfer phenomena, and the study of nuclear effects in deep inelastic 
scattering. The ability to run with both hydrogen and heavier targets will provide rich data 
sets useful for quantitative studies of nuclear models. The study of A polarization both in 
the target and in the fragmentation regions will help clarify the intriguing problem of spin 
transfer. 

Although the neutrino beam energies are well below any reasonable threshold for new 
physics, the large statistics makes it possible to search for physics beyond the Standard 
Model. The high intensity neutrino beam allows a search for the production and decay of 
neutral heavy leptons with mixing angle sensitivity two orders of magnitude better than 
present limits in the 30-80 MeV range. The exchange of new gauge bosons decoupled from 
the first generation of quarks and leptons can be seen via enhancements of the inclusive 
charm production rate, with a sensitivity well beyond the present limits. A novel neutrino 
magnetic moment search technique that uses oscillating magnetic fields at the neutrino 
beam source could discover large neutrino magnetic moments predicted by some theories. 
Rare lepton-flavor- violating decays of muons in the ring could be tagged in the deep inelastic 
scattering final states through the detection of wrong-sign electrons and muons, or of prompt 
taus. 



D. Physics that can be done with Intense Cold Muon Beams 

Experimental studies of muons at low and medium energies have had a long and distin- 
guished history, starting with the first search for muon decay to electron plus gamma-ray j3| , 
and including along the way the 1957 discovery of the nonconservation of parity, in which 
the g value and magnetic moment of the muon were first measured Q]. The years since 
then have brought great progress: limits on the standard-model-forbidden decay fi — > 
have dropped by nine orders of magnitude, and the muon anomalous magnetic moment 
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<V = (g^ — 2)/2 has yielded one of the more precise tests (~ 1 ppm) of physical theory [85]. 
The front end of a Neutrino Factory has the potential to provide ~ 10 21 muons per year, 



five orders of magnitude beyond the most intense beam currently available 163]. 

Such a facility could enable a rich variety of precision measurements. In the area of low 
energy muon physics a majority of experiments with a high physics potential is limited at 
present by statistics. The list of conceivable projects includes (see Table ITT]) : 

• precise determinations of the properties characterizing the muon, which are the mass 
m M , magnetic moment magnetic anomaly a^, charge and lifetime r M , 

• measurements the muon decay parameters (Michel parameters), 

• CPT tests from a comparison of fi~ and /i + properties, 

• measurements of fundamental constants of general importance (e.g. the electromag- 
netic fine structure constant a or the weak interaction Fermi constant Gf) 

• sensitive searches for physics beyond the Standard Model either through measuring dif- 
ferences of muon parameters from Standard Model predictions or in dedicated searches 
for rare and forbidden processes, such as fi — > e^, fi — > eee, fi~N — > e~N conversion 
and muonium-antimuonium (M — M) conversion or searches for a permanent electric 
dipole moment of the particle, 

• searches for P and T violation in muonic atoms, 

• precise determinations of nuclear properties in muonic (radioactive) atoms, 

• applications in condensed matter, thin films and at surfaces, 

• applications in life sciences, and 

• muon catalyzed fusion(/iCF). 

aeen made by a CERN study group, 



A detailed evaluation of the possibilities has recently_ 
which assumed a facility with a 4 MW proton driver 189 

In the search for "forbidden" decays, Marciano [9Cj has suggested that muon Lepton 
Flavor Violation (LFV) (especially coherent muon-to-electron conversion in the field of a 
nucleus) is the "best bet" for discovering signatures of new physics using low-energy muons. 
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The MECO experiment |9JJ proposed at BNL offers, through a novel detector concept, very 
high sensitivity and some four orders of magnitude improvement over the current limits from 
PSI [92] . At a future high muon flux facility, such as the Neutrino Factory, this could be 
improved further by 1-2 orders of magnitude. 

The search for \x — > is also of great interest. The MEGA experiment recently set 
an upper limit B{fj,+ -> e+ 7 ) < 1.2 x 10" 11 Q- 

Ways to extend sensitivity to the 10 
level have not only been discussed 9^] but also have lead to an active proposal at PSI 
The experiment aims for an improvement of three orders of magnitude over MEGA which 
had systematics limitations. The /i-to-e-conversion approach has the additional virtue of 
sensitivity to new physics that does not couple to the photon. 

An observation of a non-zero Electric Dipole Moment (EDM) of the muon, d M , could prove 
equally exciting; This has generated a Letter of Intent to observe d^, which proposes 
to use the the large electric fields associated with relativistic particles in a magnetic storage 
ring. As CP violation enters in the quark sector starting with the second generation, the 
muon is a particularly valuable probe in this regard, despite the already low limits for 
electrons. Moreover, there exist some models in which the electric dipole moment scales 
stronger than linearly 97 1. 

It is worth noting that for searches of rare muon decays and for d^ that the standard model 
predictions are orders of magnitude below the presently established limits. Any observation 
which can be shown to be not an artefact of the experimental method or due to background 
would therefore be a direct sign of new physics. 

There are three experiments going on currently to improve the muon lifetime 
Note that the Fermi coupling constant Gp is derived from a measurement of t m . The efforts 
are therefore worthwhile whenever experimental conditions allow substantial improvement. 
One should however be aware that a comparison with theory in this channel is presently 
dominated by theoretical uncertainties. 

In the case of precision measurements (r^, a^, etc.), new-physics effects appear as small 
corrections arising from the virtual exchange of new massive particles in loop diagrams. In 
contrast, LFV and EDMs are forbidden in the standard model, thus their observation at 
any level would constitute evidence for new physics. 

The current status and prospects for advances in these areas are shown in Table ILT1 which 
lists present efforts in the field and possible improvements at a Neutrino Factory or Muon 
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Type of 
Experiment 


Physics Issues 


Possible 
experiments 


Previously 
established 
accuracy 


Present activities 
(proposed accuracy) 


Projected for 
SMS @ CERN 


"Classical" 
rare & 

forbidden 
decays 


Lepton number violation; 
searches for new physics: 
SUSY, L-R Symmetry, 
R-parity violation,. . . 


fi~N -> e~N 
H e7 
—*■ eee 


6.1 ■ icr 13 

1.2 ■ 10- 11 

1.0 ■ io~ 12 

8.1 ■ 10~ u 


PSI, proposed BNL (5 ■ 10~ 17 ) 
Proposed PSI (1 ■ 10~ 14 ) 
completed 1985 PSI 
completed 1999 PSI 


< io- 18 

< ur 15 

< io- 16 

< io- 13 


Muon Decays 


Gjr; searches for new physics; 
Michel parameters 


non (V - A) 


18 ■ 10~ 6 
typ. few 10 — 3 


PSI (2x), RAL (1 ■ IO" 6 ) 
PSI, TRIUMF (1 ■ HT 3 ) 


< 10~ 7 

< io- 4 


Muon 
Moments 


Standard model tests; new 
physics; CPT tests T- resp. 
CP-violation in 2nd lepton 
generation 




1.3 ■ 10~ 6 
3.4- 10~ 19 e cm 


BNL (3.5 ■ 10~ 7 ) 
proposed BNL (10~ 24 e cm) 


< IO" 7 
< 5 ■ 10~ 26 e cm 


Muonium 
Spectroscopy 


Fundamental constants, fAu, 
WW) ot; weak interactions; 
muon charge 


M H FS 

M\ s 2s 


12 ■ 10~ 9 

1 ■ io- 9 


completed 1999 LAMPF 
completed 2000 RAL 


5 ■ IO" 9 

< io- 11 


Muonic 
Atoms 


Nuclear charge radii; 
weak interactions 


fi~ atoms 


depends 


PSI, possible CERN 
(< r p > to 10~ 3 ) 


new nuclear 
structure 


Condensed 
Matter 


surfaces, catalysis, bio 
sciences. . . 


surface /iSR 


n/a 


PSI, RAL (n/a) 


high rate 



TABLE II: Experiments which could beneficially take advantage of the intense future stopped 
muon source. The numbers were worked out for scenarios at a future Stopped Muon Source (SMS) 



of a neutrino factory at CERN 



89l | . They are based on a muon flux of 10 particles per annum 



in which beam will be available for 10 7 s. Typical beam requirements are given in Table ITTT1 



Collider facility. The beam parameters necessary for the expected improvements are listed 
in Table nm 

It is worth recalling that LFV as a manifestation of neutrino mixing is suppressed as 
(5m 2 ) 2 /my V and is thus entirely negligible. However, a variety of new-physics scenarios 
predict observable effects. Table HVl lists some examples of limits on new physics that would 
be implied by nonobservation of /i-to-e conversion (u~ N — > e~N) at the 10 -16 level |ooj | . 

The muon magnetic anomaly (muon g-2 value [10ll |) has been measured recently at the 
Brookhaven National Laboratory (BNL) with 0.7 ppm accuracy At present, no definite 
statement can be made whether this result agrees or disagrees with standard theory, which 
is sensitive to electroweak corrections. The theory has recently come under severe scrutiny 
and in particlar an error was found in the calculation of hadronic light by light scattering 



102]. The theoretical calculations are being improved upon, and with more data, there is 
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TABLE III: Beam requirements for new muon experiments. We show the needed muonic charge 
g M and the minimum of the total muon number j I^dt above which significant progress can be 
expected in the physical interpretation of the experiments. Measurements which require pulsed 
beams are sensitive to the muon suppression Io/I m between pulses of length 5T and separation 
AT. Most experiments require energies up to 4 MeV corresponding to 29 MeV/c momentum. Thin 
targ ets, respectively storage ring acceptances, demand rather small momentum spreads Ap^/p^ 



Experiment 




jl^dt 




5T 


AT 














[ns] 


[ns] 


[MeV] 


[%] 


H~N -> e~iV 


- 


10 19 


< io- 9 


< 100 


> 1000 


< 20 


1...5 


H -> e7 


+ 


10 16 


n/a 


continuous 


continuous 


1...4 


1...5 


fi — > eee 


+ 


10 15 


n/a 


continuous 


continuous 


1...4 


1...5 


fi + e~ — > u~e + 


+ 


10 16 


< io- 4 


< 1000s 


> 20000 


1...4 


1...2 




+ 


10 13 


< 1(T 4 


< 100 


> 20000 


4 


1...10 


non(V — A) 


± 


10 13 


n/a 


continuous 


continuous 


4 


1...5 




± 


10 15 


< io- 7 


< 50 


> 10 6 


3100 


10^ 4 




± 


10 16 


< 1(T 6 


< 50 


> 10 6 


<1000 


< IO" 5 


Mhfs 


+ 


10 15 


< io- 4 


< 1000 


> 20000 


4 


1...3 


M U 2s 


+ 


10 14 


< 1(T 3 


< 500 


> 10 6 


1...4 


1...2 


atoms 




10 14 


< IO" 3 


< 500 


> 20000 


1...4 


1...5 


condensed matter 


± 


10 14 


< IO" 3 


< 50 


> 20000 


1...4 


1...5 


(incl.bio sciences) 

















a goo d cha nce that this might eventually lead to evidence for beyond the standard model 



effects 



1031 ] . The final goal of the experimental precision is 0.35 ppm for the current set 



of experiments. This value could be improved by an order of magnitude at a Neutrino 
Factory, provided cold muons of energy 3.1 GeV are made available. This could further 
spur more accurate theoretical calculations that improve u pon c ontributions from hadronic 



vacuum polarization and hadronic light by light scattering 



10J]. In addition, the muon g-2 
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TABLE IV: Some examples of new physics probed by the nonobservation of \i — > e conversion at 
the 10" 16 level (from |aj|). 



New Physics 


Limit 


Heavy neutrino mixing 


KnVcn\ 2 < 10- 12 


Induced Zfie coupling 


9z, e < 10~ 8 


Induced H^ie coupling 


g Hfie < 4 x 1(T 8 


Compositeness 


A c > 3,000TeV 



experiments at CERN have provided the best test of CPT invariance at a level of 2 ■ 10 



-22 



which is more th an t hree orders of magnitude better than the mostly quoted result K° — K° 
mass difference |l0f| . A g — 2 measurement at the Neutrino Factory front end that uses 
muons of both charges would lead to further improvement in these CPT limits. 

Precision studies of atomic electrons have provided notable tests of QED ( e.g, the Lamb 
shift in hydrogen) and could in principle be used to search for new physics were it not for 
nuclear corrections. Studies of muonium (/i + e~) are free of such corrections since it is a 
purely leptonic system. Muonic atoms can also yield new information complementary to 
that obtained from electronic atom s. A number of p ossibilities have been enumerated by 
Kawall et al. 106], Jungmann jlO t| and Molzon |l08| . 



By making measurements on the muonium system, for instance, one can produce precise 
measurements of the fundamental constants and also do sensitive searches fo r new physics. 
The muonium ground state hyperfine structure has been measured to 12 ppb jiod ] and cur- 
rently furnishes the most sensitive test of the relativistic two-body bound state in QED 107| . 
The precision could be further improved significantly with increased statistics. The theoret- 
ical error is 120 ppb. The uncertainty arising from the muon mass is five times larger than 
that from calculations. If one assumes the theory to be correct, the muon-electron mass 
ratio can be extracted to 27 ppb. A precise value for the electromagnetic fine structure con- 
stant a can be extracted. Its good agreement with the number extracted from the electron 
magnetic anomaly must be viewed as the best test of the internal consistency of QED, as 
one case involves bound state QED and the other that of free particles. The Zeeman effect of 
the muonium hyperfine structure allows the best direct measurement of the muon magnetic 
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moment, respectively its mass, to 120 ppb, improved by higher-precision measurements in 
muonium and muon spin resonance. These are also areas in which the Neutrino Fac tory 
front end could contribute. Laser spectroscopy of the muonium ls-2s transition llOj ] has 
resulted in a precise value of the muon mass as well as the testing of the muon-electron 
charge ratio to about 2 • 1CT 9 . This is by far the best test of charge equality in the first two 
generations. 

The search for muonium-antimuonium conversion had been prop osed by Pontecorvo 
three years before the systemwas first produced by Hughes et al. Several new- 

physics models allow viol ation of lepton family number by two units. The current limit 
is R g = Gc/Gp < 0.0030 |l!2| . where Gc is the new-physics coupling constant, and Gp the 
Fermi coupling constant. This sets a lower limit of 2.6TeV/c 2 (90% C.L.) on the mass of a 
grand-unified dileptonic gauge bos on a nd also strongly disfavours models with heavy lepton 
seeded radiative mass generation 1121. The search for muonium-antimuoniu m con version 

n 

has by far the strongest gain in sensitivity of all rare muon decay experiments 107J . 

The high intensity proton machine needed for the Neutrino Factory can also find use as a 
new generation isotope facility which would have much higher rates compared to the present 
ISOLDE facility at CERN. Nucleids yet not studied could be produced at quantities which 
allow precision investigations of their properties js^]. The measurements of muonic spectra 
can yield most precise values for the charge radii of nuclei as well as other ground state 
properties such as moments and even B(E2) transition strengths for even-even nuclei. An 
improved understanding of nuclear structure can be expected which may be of significance for 
interpreting various neutrino experiments, rare decays involving nuclei, and nuclear capture. 
An urgent need exists for accurate charge and neutron radii of Francium and Radium isotopes 
which are of interest for atomic parity violation research and EDM searches in atoms and 
nuclei. 

Muonic x-ray experiments generally promise higher accuracy for most of these quantities 
compared to electron scattering, particularly because the precision of electron scattering 
data depends on the location of the minimum of the cross section where rates are naturally 
low. In principle, for chains of isotopes charge radii can be inferred from isotope shift 
measurements with laser spectroscopy. However, this gives only relative information. For 
absolute values, calibration is necessary and has been obtained in the past for stable nuclei 
from muonic spectra. In general, two not too distant nuclei are needed for a good calibration. 
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The envisa ged experimental approaches include i) the technique pioneered by Nagamine 
and Strasser 113J, which involves cold films for keeping radioactive atoms and as a host 
material in which muon transfer takes place; ii) merging beams if radioactive ions and of 
muons; and iii) trapping of exotic isotopes in a Penning trap which is combined with a 
cyclo tron trap. Large formation rates can be expected from a setup conta ining a Penning 



trap 



114j . the magnetic field of which serves also as a cyclotron muon trap 1151 ] . For muon 



energies in the range of electron binding energies the muon capture cross sections grow to 
atomic values, efficient atom production results at the rate of approxi mate ly 50 Hz. It 
should be noted that antiprotonic atoms could be produced similarly |llfi| ] and promise 
measurements of neutron distributions in nuclei. 



E. Physics potential of a Low energy Muon Collider operating as a Higgs Factory 

Muon colliders 



117 . Ill8l l have a number of unique features that make them attractive 

n 

candidates for future accelerators [9J. The most important and fundamental of these derive 
from the large mass of the muon in comparison to that of the electron. The synchrotron 
radiation loss in a circular accelerator goes as the inverse fourth power of the mass and 
is two billion times less for a muon than for an electron. Direct s channel coupling to 
the higgs boson goes as the mass squared and is 40,000 greater for the muon than for 
the electron. This leads to: a) the possibility of extremely narrow beam energy spreads, 
especially at beam energies below 100 GeV; b) the possibility of accelerators with very high 
energy; c) the possiblity of employing storage rings at high energy; d) the possibility of using 
decays of accelerated muons to provide a high luminosity source of neutrinos as discussed 
in Section 111 A 41 e) increased potential for probing physics in which couplings increase with 
mass (as does the SM hsMff coupling) . 

The relatively large mass of the muon compared to the mass of the electron means that 
the coupling of Higgs bosons to fi + n~ is very much larger than to e + e~, implying much 
larger s-channel Higgs production rates at a muon collider as compared to an electron 
collider. For Higgs bosons with a very small (MeV-scale) width, such as a light SM Higgs 
boson, production rates in the s-channel are further enhanced by the muon collider's ability 
to achieve beam energy spreads comparable to the tiny Higgs width. In addition, there is 
little beamstrahlung, and the beam energy can be tuned to one part in a million through 
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continuous spin-rotation measurements 



119l | . Due to these important qualitative differences 



between the two types of machines, only muon colliders can be advocated as potential s- 
channel Higgs factories capable of determining the mass and decay width of a Higgs boson 



to very high precision 



120. 



12l( | . High rates of Higgs production at e + e colliders rely on 



substantial VV Higgs coupling for the Z+Higgs (Higgstrahlung) or WW ^Higgs (WW 
fusion) reactions. In contrast, a fi + fi~ collider can provide a factory for producing a Higgs 
boson with little or no VV coupling so long as it has SM-like (or enhanced) couplings. 

Of course, there is a tradeoff between small beam energy spread, SE/E = R, and lumi- 
nosity. Current estimates for yearly integrated luminosities (using C = 1 x 10 32 cm -2 s _1 
as implying L = 1 fb _1 /yr) are: L ycar > 0.1,0.22,1 fb _1 at y/s ~ 100 GeV for beam 
energy resolutions of R = 0.003%, 0.01%, 0.1%, respectively; L ycaT ~ 2,6,10 fb _1 at 
s/s ~ 200,350,400 GeV, respectively, for R ~ 0.1%. Despite this, studies show that for 
small Higgs width the s-channel production rate (and statistical significance over back- 
ground) is maximized by choosing R to be such that < r^ ot . In particular, in the SM 
context for mh SM ~ 110 GeV this corresponds to R ~ 0.003%. 



If the rrih ~ 115 GeV LEP signal is real, or if the interpretation of the precision elec- 
troweak data as an indication of a light Higgs boson (with substantial VV coupling) is valid, 
then both e + e~ and colliders will be valuable. In this scenario the Higgs boson would 

have been discovered at a previous higher energy collider (even possibly a muon collider run- 
ning at high energy), and then the Higgs factory would be built with a center-of-mass energy 
precisely tuned to the Higgs boson mass. The most likely scenario is that the Higgs boson 
is discovered at the LHC via gluon fusion (gg — > H) or perhaps earlier at the Tevatron via 
associated production (qq — > WH, ttH), and its mass is determined to an accuracy of about 
100 MeV. If a linear collider has also observed the Higgs via the Higgs-strahlung process 
(e + e~ — > ZH), one might know the Higgs boson mass to better than 50 MeV with an inte- 
grated luminosity of 500 fb _1 . The muon collider would be optimized to run at y/s ~ m#, 
and this center-of-mass energy would be varied over a narrow range so as to scan over the 
Higgs resonance (see Fig. ITUl below). 
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1. Higgs Production 



The production of a Higgs boson (genetic ally denot ed h) in the s-channel with interesting 



rates is a unique feature of a muon collider 



120 



12 1| . The resonance cross section is 



,2 • 



(32) 



[s-mlf + ml (ry ; 

In practice, however, there is a Gaussian spread (cx^) to the center-of-mass energy and one 
must compute the effective s-channel Higgs cross section after convolution assuming some 
given central value of \fs: 
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It is convenient to express in terms of the root-mean-square (rms) Gaussian spread of 
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FIG. 10: Number of events and statistical errors in the bb final state as a function of y/s in the 
vicinity of rrih SAI = 110 GeV, assuming R = 0.003%, and eL = 0.00125 fb _1 at each data point. 
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the energy of an individual beam, R: 

^ = < 2 MeV > (^) (lWllv) ' < 35 » 

From Eq. (|32|) . it is apparent that a resolution cr^ < r}°* is needed to be sensitive to the Higgs 
width. Further, Eq. (|3*4*j) implies that oc 1/ct^ for cr^ > r^ ot and that large event rates 
are only possible if is not so large that BF(h — * /i/2) is extremely suppressed. The width 
of a light SM-like Higgs is very small ( e.g, a few MeV for m/ lSM ~ 110 GeV), implying the 
need for R values as small as ~ 0.003% for studying a light SM-like h. Figure ITU1 illustrates 
the result for the SM Higgs boson of an initial centering scan over y/s values in the vicinity 
of mh SM = 110 GeV. This figure dramatizes: a) that the beam energy spread must be very 
small because of the very small rj*^ (when m/ lgM is small enough that the WW* decay 
mode is highly suppressed); b) that we require the very accurate in situ determination of the 
beam energy to one part in a million through the spin precession of the muon noted earlier 
in order to perform the scan and then center on = rrih SM with a high degree of stability. 
If the h has SM-like couplings to WW, its width will grow rapidly for > 2m>w and its 
s-channel production cross section will be severely suppressed by the resulting decrease of 
BF(h — > fifi). More generally, any h with SM-like or larger hfifi coupling will retain a large 
s-channel production rate when > 2myj/ only if the hWW coupling becomes strongly 
suppressed relative to the IismWW coupling. 

The general theoretical prediction within supersymmetric models is that the lightest 
supersymmetric Higgs boson h° will be very similar to the hsu when the other Higgs bosons 
are heavy. This 'decoupling limit' is very likely to arise if the masses of the supersymmetric 
particles are large (since the Higgs masses and the superparticle masses are typically similar 
in size for most boundary condition choices). Thus, h° rates will be very similar to hsM rates. 
In contrast, the heavier Higgs bosons in a typical supersymmetric model decouple from VV 
at large mass and remain reasonably narrow. As a result, their s-channel production rates 
remain large. 

For a SM-like h, at yfs = w 115 GeV and R = 0.003%, the bb rates are 

signal w 10 4 events x L(fb _1 ) , (36) 
background w 10 4 events x L(fb _1 ) . (37) 
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2. What the Muon Collider Adds to LHC and LC Data 

An assessment of the need for a Higgs factory requires that one detail the unique capa- 
bilities of a muon collider versus the other possible future accelerators as well as comparing 
the abilities of all the machines to measure the same Higgs properties. Muon colliders, and a 
Higgs factory in particular, would only become operational after the LHC physics program 
is well-developed and, quite possibly, after a linear collider program is mature as well. So 
one important question is the following: if a SM-like Higgs boson and, possibly, important 
physics beyond the Standard Model have been discovered at the LHC and perhaps studied 
at a linear collider, what new information could a Higgs factory provide? The s-channel 
production process allows one to determine the mass, total width, and the cross sections 
Whi^ \i~ — > h — > X) for several final states X to very high precision. The Higgs mass, total 
width and the cross sections can be used to constrain the parameters of the Higgs sector. 
For example, in the MSSM their precise values will constrain the Higgs sector parameters 
mpfi and tan (3 (where tan (3 is the ratio of the two vacuum expectation values (vevs) of the 
two Higgs doublets of the MSSM). The main question is whether these constraints will be a 
valuable addition to LHC and LC constraints. The expectations for the luminosity available 
at linear colliders has risen steadily. The most recent studies assume an integrated lumi- 
nosity of some 500 fb _1 corresponding to 1-2 years of running at a fewxlOO fb _1 per year. 
This luminosity results in the production of greater than 10 4 Higgs bosons per year through 
the Bjorken Higgs-strahlung process, e + e~ — *> Zh, provided the Higgs boson is kinematically 
accessible. This is comparable or even better than can be achieved with the current machine 
parameters for a muon collider operating at the Higgs resonance; in fact, recent studies have 
described high-luminosity linear colliders as "Higgs factories," though for the purposes of 
this report, we will reserve this term for muon colliders operating at the s-channel Higgs 
resonance. A linear collider with such high luminosity can certainly perform quite accu- 
rate measurements of certain Higgs par amet ers, such as the Higgs mass, couplings to gauge 
bosons and couplings to heavy quarks |l22j |. Precise measurements of the couplings of the 
Higgs boson to the Standard Model particles is an important test of the mass generation 
mechanism. In the Standard Model with one Higgs doublet, this coupling is proportional 
to the particle mass. In the more general case there can be mixing angles present in the 
couplings. Precision measurements of the couplings can distinguish the Standard Model 
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Higgs boson from that from a more general model and can constrain the parameters of a 
more general Higgs sector. 

TABLE V: Achievable relative uncertainties for a SM-like rrih = 110 GeV for measuring the Higgs 
boson mass and total width for the LHC, LC (500 fbr 1 ), and the muon collider (0.2 fb 1 ). 





LHC 


LC n+p- 




9 x 10~ 4 


3 x 10" 4 1 - 3 x 10~ 6 


ptot 

1 h 


> 0.3 


0.17 0.2 



The accuracies possible at different colliders for measuring rrih and T^ * of a SM-like h 
with rrih ~ 110 GeV are given in Table IVl Once the mass is determined to about 1 MeV at 



the LHC and/or LC, the muon collider would employ a three-point fine scan |l20j near the 
resonance peak. Since all the couplings of the Standard Model are known, is known. 

Therefore a precise determination of T h ot is an important test of the Standard Model, and 
any deviation would be evidence for a nonstandard Higgs sector. For a SM Higgs boson with 
a mass sufficiently below the WW* threshold, the Higgs total width is very small (of order 
several MeV), and the only process where it can be measured directly is in the s-channel at a 
muon collider. Indirect determinations at the LC can have higher accuracy once rrih is large 
enough that the WW* mode rates can be accurately measured, requiring rrih > 120 GeV. 
This is because at the LC the total width must be determined indirectly by measuring a 
partial width and a branching fraction, and then computing the total width, 

T(h^X) 



tot 



BR(h -> X) ' 



(38) 



for some final state X. For a Higgs boson so light that the WW* decay m ode is not useful, 
the total width measurement would probably require use of the 77 decays jl23| . This would 
require information from a photon collider as well as the LC and a small error is not possible. 
The muon collider can measure the total width of the Higgs boson directly, a very valuable 
input for precision tests of the Higgs sector. 

To summarize, if a Higgs is discovered at the LHC or possibly earlier at the Fermilab 
Tevatron, attention will turn to determining whether this Higgs has the properties expected 
of the Standard Model Higgs. If the Higgs is discovered at the LHC, it is quite possible that 
super symmetric states will be discovered concurrently. The next goal for a linear collider or a 

44 



muon collider will be to better measure the Higgs boson properties to determine if everything 
is consistent within a supersymmetric framework or consistent with the Standard Model. 
A Higgs factory of even modest luminosity can provide uniquely powerful constraints on 
the parameter space of the supersymmetric model via its very precise measurement of the 
light Higgs mass, the highly accurate determination of the total rate for — > h° — ► bb 
(which has almost zero theoretical systematic uncertainty due to its insensitivity to the 
unknown m;, value) and the moderately accurate determination of the h 0, s total width. In 
addition, by combining muon collider data with LC data, a completely model-independent 
and very precise determination of the h°n + n~ coupling is possible. This will provide another 
strong discriminator between the SM and the MSSM. Further, the h°fi + fi~ coupling can be 
compared to the muon collider and LC determinations of the h°T + r~ coupling for a precision 
test of the expected universality of the fermion mass generation mechanism. 



F. Physics Potential of a High Energy Muon Collider 



Once one learns to cool muons, it becomes possible to build muon cp 
3 TeV in the center of mass that fit on an existing l aboratory site 
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liders with energies of 
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ener gies, it becomes possible to measure the W mass 

Il2 t| with high accuracy, by scanning the thresholds of these particles and making use of the 
excellent energy resolution of the beams. We consider further here the ability of a higher 
energy muon collider to scan higher-lying Higgs like objects such as the H° and the A in 
the MSSM that may be degenerate with each other. 



1. Heavy Higgs Bosons 

As discussed in the previous section, precision measurements of the light Higgs boson 
properties might make it possible to not only distinguish a supersymmetric boson from a 
Standard Model one, but also pinpoint a range of allowed masses for the heavier Higgs 
bosons. This becomes more difficult in the decoupling limit where the differences between 
a supersymmetric and Standard Model Higgs are smaller. Nevertheless with sufficiently 
precise measurements of the Higgs branching fractions, it is possible that the heavy Higgs 
boson masses can be inferred. A muon collider light-Higgs factory might be essential in this 
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process. In the context of the MSSM, m^o can probably be restricted to within 50 GeV or 
better if m^o < 500 GeV. This includes the 250 — 500 GeV range of heavy Higgs boson 
masses for which discovery is not possible via H°A° pair production at a y/s = 500 GeV 
LC. Further, the A and H° cannot be detected in this mass range at either the LHC or LC 
in bbH°,bbA production for a wedge of moderate tan/5 values. (For large enough values 
of tan/3 the heavy Higgs bosons are expected to be observable in bbA°,bbH° production at 
the LHC via their r + r~ decays and also at the LC.) A muon collider can fill some, perhaps 
all of this moderate tan/9 wedge. If tan/3 is large, the [i + fi~H° and fi + fi~A° couplings 
(proportional to tan/3 times a SM-like value) are enhanced, thereby leading to enhanced 
production rates in collisions. The most efficient procedure is to operate the muon 

collider at maximum energy and produce the H° and A (often as overlapping resonances) 
via the radiative return mechanism. By looking for a peak in the bb final state, the H° 
and A can be discovered and, once discovered, the machine can be set to m^o or m#o 
and factory-like precision studies pursued. Note that the A and H° are typically broad 
enough that R = 0.1% would be adequate to maximize their s-channel production rates. In 
particular, T ~ 30 MeV if the it decay channel is not open, and T ~ 3 GeV if it is. Since 
R = 0.1% is sufficient, much higher luminosity (L ~ 2 — 10 fb~ /yr) would be possible as 
compared to that for R = 0.01% — 0.003% required for studying the h°. 

In short, for these moderate tan/3-m^o > 250 GeV scenarios that are particularly difficult 
for both the LHC and the LC, the muon collider would be the only place that these extra 
Higgs bosons can be discovered and their properties measured very precisely. 

In the MSSM, the heavy Higgs bosons are largely degenerate, especially in the decoupling 
limit where they are heavy. Large values of tan /3 heighten this degeneracy. A muon collider 
with sufficient energy resolution might be the only possible means for separating out these 
states. Examples showing the H and A resonances for tan/3 = 5 and 10 are shown in Fig.fTTl 
For the larger value of tan (3 the resonances are clearly overlapping. For the better energy 
resolution of R = 0.01%, the two distinct resonance peaks are still visible, but become 
smeared out for R = 0.06%. 

Once muon colliders of these intermediate energies can be built , higher energies such as 
3-4 TeV in the center of mass become feasible. Muon colliders with these energies will be 
complementary to hadron colliders of the SSC class and above. The backgroun d ra diation 



from neutrinos from the muon decay becomes a problem at ~ 3 TeV in the CoM 
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FIG. 11: Separation of A and H signals for tan/3 = 5 and 10. From Ref. 12C] 



for ameliorating this problem have been discussed and include optical stochastic cooling to 
reduce the number of muons needed for a given luminosity, elimination of straight sections 
via wigglers or undulators, or special sites for the collider such that the neutrinos break 
ground in uninhabited areas. 

III. NEUTRINO FACTORY 



In this Section we describe the various components of a Neutrino Factory, based on 
the most recent Feasibility Study (Study II) Q that was carried out jointly by BNL and 
the MC. We also describe the stages that could be constructed incrementally to provide 
a productive physics program that evolves eventually into a full-fledged Neutrino Factory. 
Details of the design described here are based on the specific scenario of sending a neutrino 
beam from Brookhaven to a detector in Carlsbad, New Mexico. More generally, however, the 
design exemplifies a Neutrino Factory for which our two Feasibility Studies demonstrated 
technical feasibility (provided the challenging component specifications are met), established 
a cost baseline, and established the expected range of physics performance. As noted earlier, 
this design typifies a Neutrino Factory that could fit comfortably on the site of an existing 
laboratory, such as BNL or FNAL. 

A list of the main ingredients of a Neutrino Factory is given below: 
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• Proton Driver: Provides 1-4 MW of protons on target from an upgraded AGS; a 
new booster at Fermilab would perform equivalently. 

• Target and Capture: A high-power target immersed in a 20 T superconducting 
solenoidal field to capture pions produced in proton-nucleus interactions. 

• Decay and Phase Rotation: Three induction linacs, with internal superconducting 
solenoidal focusing to contain the muons from pion decays, that provide nearly non- 
distorting phase rotation; a "mini-cooling" absorber section is included after the first 
induction linac to reduce the beam emittance and lower the beam energy to match 
the downstream cooling channel acceptance. 

• Bunching and Cooling: A solenoidal focusing channel, with high-gradient rf cavities 
and liquid hydrogen absorbers, that bunches the 250 MeV/c muons into 201.25 MHz 
rf buckets and cools their transverse normalized rms emittance from 12 mm-rad to 2.7 
mm-rad. 

• Acceleration: A superconducting linac with solenoidal focusing to raise the muon 
beam energy to 2.48 GeV, followed by a four-pass superconducting RLA to provide a 
20 GeV muon beam; a second RLA could optionally be added to reach 50 GeV, if the 
physics requires this. 

• Storage Ring: A compact racetrack-shaped superconducting storage ring in which 
~35% of the stored muons decay toward a detector located about 3000 km from the 
ring. 

A. Proton Driver 

The proton driver considered in Study II is an upgrade of the BNL Alternating Gradient 
Synchrotron (AGS) and uses most of the existing components and facilities; parameters are 
listed in Table IVT1 To serve as the proton driver for a Neutrino Factory, the existing booster 
is replaced by a 1.2 GeV superconducting proton linac. The modified layout is shown in 
Fig. El The AGS repetition rate is increased from 0.5 Hz to 2.5 Hz by adding power 
supplies to permit ramping the ring more quickly. No new technology is required for this — 
the existing supplies are replicated and the magnets are split into six sectors rather than the 
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To RHIC 




FIG. 12: (Color)AGS proton driver layout. 

two used presently. The total proton charge (10 14 ppp in six bunches) is only 40% higher 
than the current performance of the AGS. However, due to the required short bunches, there 
is a large increase in peak current and concomitant need for an improved vacuum chamber; 
this is included in the upgrade. The six bunches are extracted separately, spaced by 20 ms, 
so that the target, induction linacs, and rf systems that follow need only deal with single 
bunches at an instantaneous repetition rate of 50 Hz (average rate of 15 Hz). The average 
proton beam power is 1 MW. A possible future upgrade to 2 xlO 14 ppp and 5 Hz could 
give an average beam power of 4 MW. At this higher intensity, a superconducting bunch 
compressor ring would be needed to maintain the rms bunch length at 3 ns. 

If the facility were built at Fermilab, th e pr oton driver would be a newly constructed 
16 GeV rapid cycling booster synchrotron |l29j . The planned facility layout is shown in 
Fig. El The initial beam power would be 1.2 MW, and a future upgrade to 4 MW would be 
possible. The Fermilab design parameters are included in Table IVII A less ambitiou s and 
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more cost-effective 8 GeV proton driver option has also been considered for Fermilab 
this too might be the basis for a proton driver design. 

B. Target and Capture 

A mercury jet target is chosen to give a high yield of pions per MW of incident proton 
power. The 1 cm diameter jet is continuous, and is tilted with respect to the beam axis. The 
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FIG. 13: (Color)FNAL proton driver layout from Ref. 



target layout is shown in Fig. El We assume that the thermal shock from the interacting 
proton bunch fully disperses the mercury, so the jet must have a velocity of 20-30 m/s 
to be replaced before the next bunch. Calculations of pion yields that reflect the det ailed 
magnetic geometry of the target area have been performed with the MARS code jisj ]. To 
avoid mechanical fatigue problems, a mercury pool serves as the beam dump. This pool 
is part of the overall target — its mercury is circulated through the mercury jet nozzle after 
passing through a heat exchanger. 

Pions emerging from the target are captured and focused down the decay channel by a 
solenoidal field that is 20 T at the target center, and tapers down, over 18 m, to a periodic 
(0.5 m) superconducting solenoid channel (B z = 1.25 T) that continues through the phase 
rotation to the start of bunching. Note that the longitudinal direction of the fields in this 
channel do not change sign from cell to cell as they do in the cooling channel. The 20 T 
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TABLE VI: Proton driver parameters for BNL and FNAL designs. 





BNL 


FNAL 


Total beam power (MW) 


1 


1.2 


Beam energy (GeV) 


24 


16 


Average beam current (//A) 


42 


72 


Cycle time (ms) 


400 


67 


Number of protons per fill 


1 x 10 14 


3 x 10 13 


Average circulating current (A) 


6 


2 


No. of bunches per fill 


6 


18 


No. of protons per bunch 


1.7 x 10 13 


1.7 x 10 12 


Time between extracted bunches (ms) 


20 


0.13 


Bunch length at extraction, rms (ns) 


3 


1 



solenoid, with a resistive magnet insert and superconducting outer coil, is similar in character 
• ;o th e higher field (up to 45 T), but smaller bore, magnets existing at several laboratories 
133| . The magnet insert is mad e wit h hollow copper conductor having ceramic insulation 
to withstand radiation. MARS simulations of radiation levels show that, with the 

shielding provided, both the copper and superconducting magnets could have a lifetime 
greater than 15 years at 1 MW. 

In Study I, the target was a solid carbon rod. At high beam energies, this implementation 
has a lower pion yield than the mercury jet, and is expected to be more limited in its ability 
to handle the proton beam power, but should simplify the target handling issues that must 
be dealt with. At lower beam energies, say 6 GeV, the yield difference between C and Hg 
essentially disappears, so a carbo n target would be a competitive option with a lower energy 
driver. Present indications |l34| are that a carbon-carbon composite target can be tailored 
to tolerate even a 4 MW proton beam power — a very encouraging result. Other alternative 
approaches, including a rotating Inconel band target, and a granular Ta target are also 
under consideration, as discussed in Study II j^j . Clearly there are several target options 
that could be used for the initial facility. 
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FIG. 14: (Color)Target, capture solenoids and mercury containment. 
C. Phase Rotation 

The function of the phase rotation section in a neutrino factory is to reduce the energy 
spread of the collected muon beam to a manageable level, allowing reasonable throughput 
in the subsequent system components. The following description refers specifically to the 
properties of the U.S. Feasibility Study 2 for a neutrino factory. The initial pions are 
produced in the mercury target with a very wide range of momenta. The momentum 
spectrum peaks around 250 MeV/c, but there is a tail of high energy pions that extends 
well beyond 1 GeV. The pions are spread in time over about 3 ns, given by the pulse 
duration of the proton driver. After the 18 m long tapered collection solenoid and an 18 m 
long drift section, where the beam is focused by 1.25 T solenoids, most of the low energy 
pions have decayed into muons. At this point the muon energy spectrum also extends over 
an approximately 1 GeV range and the time spectrum extends over approximately 50 ns. 
However, there is a strong correlation between the muon energy and time that can be used 
for "phase rotation" . 

In the phase rotation process an electric field is applied at appropriate times to decelerate 
the leading high energy muons and to accelerate the trailing low energy ones. Since the 
bunch train required by a neutrino factory can be very long, it is possible to minimize the 
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TABLE VII: Properties of the induction linacs used in Feasibility Study 2. 
Induction Linac 12 3 



Length m 100 80 100 

Peak gradient MV/m 1.5 -1.5 1.0 
Pulse FWHM ns 250 100 380 
Pulse start offset ns 55 55 



energy spread using induction linacs. The induction linac consists of a simple non-resonant 
structure, where the drive voltage is applied to an axially symmetric gap that encloses a 
toroidal ferromagnetic material. The change in flux in the magnetic core induces an axial 
electric field that provides particle acceleration. The induction linac is typically a low 
gradient structure that can provide acceleration fields of varying shapes and time durations 
from tens of nanoseconds to several microseconds. Some properties of the induction linacs 
are given in Table IVII1 
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FIG. 15: (Color)Cross section of the induction cell and transport solenoids. 



Three induction linacs are used in a system that reduces distortion in the phase-rotated 
bunch, and permits all induction units to operate with u nipol ar pulses. The induction units 



are similar to those being built for the DARHT project 



135|| . The 1.25 T beam transport 
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solenoids are placed inside the induction cores in order to avoid saturating the core material, 
as shown in Fig. [T5l 

Between the first and second induction linacs two liquid hydrogen absorbers (each 1.7 m 
long and 30 cm radius) are used to (1) provide some initial cooling of the transverse emittance 
of the muon beam and (2) lower the average momentum of the beam to match better the 
downstream cooling channel acceptance. This process is referred to as "mini-cooling". The 
direction of the solenoid magnetic field is reversed between the two absorbers. The presence 
of material in the beam path destroys the conservation of canonical angular momentum 
that occurs when a particle enters and leaves a solenoid in vacuum. The build-up of this 
angular momentum would eventually lead to emittance growth. However, this growth can 
be minimized by periodically reversing the direction of the field. 

The beam at the end of the phase rotation section has an average momentum of about 
250 MeV/c and an rms fractional energy spread of ps4.4%. Figure [HD shows the evolution 
of the beam distribution in the phase rotation section. 

D. Buncher 

The long beam pulse (400 ns) after the phase rotation is then bunched at 201.25 MHz prior 
to cooling and acceleration at that frequency. The bunching is done in a lattice identical 
to that at the start of the cooling channel, and is preceded by a matching section from 
the 1.25 T solenoids into this lattice. The bunching has three stages, each consisting of rf 
(with increasing acceleration) followed by drifts (with decreasing length). In the first two 
rf sections, second-harmonic 402.5 MHz rf is used together with the 201.25 MHz primary 
frequency to improve the capture efficiency. The 402.5 MHz cavities are designed to fit into 
the bore of the focusing solenoids, in the location corresponding to that of the liquid hydrogen 
absorber in the downstream cooling channel. Their aperture radius for the 402.5 MHz 
cavities is 20 cm at the IRIS, while that of the 201.25 MHz cavities is 25 cm. The gradients 
on axis in the cavities are 6.4 MV/m for the 402.5 MHz cavities, and range from 6 to 8 MV/m 
for the 201.25 MHz cavities. The resulting bunches fill the 201.25 MHz stationary rf bucket. 
Figure El shows the evolution of the longitudinal distribution in the buncher. 
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FIG. 16: Evolution of the beam distribution in the phase rotation section. The graphs show the 
distribution before the phase rotation, after the first induction linac (top row, left to right), after 
mini-cooling, and after the second and third induction linacs (bottom row). 

E. Cooling 

The transverse emittance of the muon beam after phase rotation and bunching must 
be reduced in order to fit into the downstream accelerators and storage ring. Ionization 
cooling is currently the only feasible option for cooling the beam within the muon lifetime. 
In ionization cooling the transverse and longitudinal momenta are lowered in the absorbers, 
but only the longitudinal momentum is restored by the rf. The following description refers 
specifically to the properties of the U.S. Feasibility Study 2 for a neutrino factory. Transverse 
emittance cooling is achieved using cooling cells that (1) lower the beam energy by 7-12 MeV 
in liquid hydrogen absorbers, (2) use 201 MHz rf cavities to restore the lost energy, and (3) 
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FIG. 17: (Color)Evolution of beam in buncher. Plots are at the beginning of the buncher (top 
left), and at the ends of the three bunching stages (top right, bottom left, and bottom right, in 
that order). 

use 3-5 T solenoids to strongly focus the beam at the absorbers. At the end of the cooling 
channel the rms normalized transverse emittance is reduced to about 2.5 mm rad. 

Each cell of the lattice contains three solenoids. The direction of the solenoidal field 
reverses in alternate cells in order to prevent the build-up of canonical angular momentum, 
as mentioned above in the discussion of mini-cooling. In analogy with the FODO lattice 
this focusing arrangement is referred to as a (S)FOFO lattice. Multiple Coulomb scattering 
together with the focusing strength determine the asymptotic limit on the transverse emit- 
tance that the cooling channel can reach. The focusing strength in the channel is tapered 
so that the angular spread of the beam at the absorber locations remains large compared to 
the characteristic spread from scattering. This is achieved by keeping the focusing strength 
inversely proportional to the emittance, i.e., increasing it as the emittance is reduced. The 
solenoidal field profile was chosen to maximize the momentum acceptance (±22%) through 
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FIG. 18: (Color)Two cells of the 1.65 m cooling lattice. 

the channel. To maintain the tapering of the focusing it was eventually necessary to reduce 
the cell length from 2.75 m in the initial portion of the channel to 1.65 m in the final portion. 
A layout of the shorter cooling cells is shown in Fig. 
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FIG. 19: The transverse (filled circles, in mm) and longitudinal (open circles, in cm) emittances, 
as a function of the distance down the cooling channel. 



Figure EH shows a simulation of cooling in this channel. The transverse emittance de- 
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FIG. 20: Muons per incident proton in the cooling channel that would fall within a normalized 
transverse acceptance of 15 mm (open circles) or 9.75 mm (filled circles). 



creases steadily along the length of the channel. This type of channel only cools transversely, 
so the longitudinal emittance increases until the rf bucket is full and then remains fairly con- 
stant as particles are lost from the bucket. A useful figure of merit for cooling at a neutrino 
factory is the increase in the number of muons that fit within the acceptance of the down- 
stream accelerators. This is shown in Fig. I^UJ At each axial position the number of muons 
is shown that fall within two acceptances appropriate to a downstream accelerator. Both 
acceptances require the muon longitudinal phase space be less than 150 mm. The density 
of particles within a normalized transverse acceptance, for example, steadily increases by 
a factor of about 3 over the channel length, clearly showing the results of cooling. The 
saturation of the yield determined the chosen channel length of 108 m. 



F. Acceleration 

The layout of the acceleration system is shown in Fig. |^ and its parameters are listed 
in Table IVIIII The acceleration system consists of a preaccelerator linac followed by a 
four-pass recirculating linac. The recirculating linac allows a reduction in the amount of rf 
required for acceleration by passing the beam through linacs multiple times. The linacs are 
connected by arcs, and a separate are is used for each pass. At low energies, however, the 
large emittance of the beam would require a much shorter cell length and larger aperture 
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FIG. 21: (Color)Accelerating system layout. 



than is desirable and needed at higher energies. This, combined with difficulties in injecting 
the large emittance and energy spread beam into the recirculating accelerator, and the loss 
of efficiency due to the phase slip at low energies, lead to the necessity for a linac that 
precedes the recirculating linac. 

A 20 m SFOFO matching section, using normal conducting rf systems, matches the beam 
optics to the requirements of a 2.87 GeV superconducting rf linac with solenoidal focusing. 
The linac is in three parts. The first part has a single 2-cell rf cavity unit per period. The 
second part, as a longer period becomes possible, has two 2-cell cavity units per period. The 
last section, with still longer period, accommodates four 2-cell rf cavity units per period. 
See Tables IIXI and |X] for details of the rf cryostructures and cavities. Figure E21 shows the 
three cryomodule types that make up the linac. 




FIG. 22: (Color)Layouts of short (top), intermediate (middle) and long (bottom) cryomodules. 
Blue lines are the SC walls of the cavities. Solenoid coils are indicated in red. 

This linac is followed by a single four-pass recirculating linear accelerator (RLA) that 
raises the energy from 2.5 GeV to 20 GeV. The RLA uses the same layout of four 2- 
cell superconducting rf cavity structures as the long cryomodules in the linac, but utilizes 
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TABLE VIII: Main parameters of the muon accelerator. 



Injection momentum (MeV/c) /Kinetic energy (MeV) 


210/129.4 


Final energy (GeV) 


20 


Initial normalized transverse acceptance (mm rad) 


15 


rms normalized transverse emittance (mm rad) 


2.4 


Initial longitudinal acceptance, ApLb/m^ (mm) 


170 


momentum spread, Ap/p 


±0.21 


bunch length, Lj, (mm) 


±407 


rms energy spread 


0.084 


rms bunch length (mm) 


163 


Number of bunches per pulse 


67 


Number of particles per bunch/per pulse 


4.4 x 10 10 /3 x 10 12 


Bunch frequency/accelerating frequency (MHz) 


201.25/ 201.25 


Average beam power (kW) 


150 



TABLE IX: Parameters for three types of linac cryomodules. 



Cavities per period 


1 


2 


4 


Period length (m) 


5 


8 


13 


Number of periods 


11 


16 


19 


Cavity type 


A 


A 


B 


Solenoid full aperture (cm) 


46 


46 


36 


Solenoid length (m) 


1 


1 


1.5 


Maximum solenoid field (T) 


2.1 


2.1 


4.2 



quadrupole triplet focusing, as indicated in Fig. [2S1 

The arcs have an average radius of 62 m, and are all in the same horizontal plane. They 
also utilize triplet focusing. There are around 120 arc cells, with 2.15 m dipoles, and triplet 
quadrupoles which are very similar to those in the linacs. The required full quadrupole 
apertures vary from 20 cm to 12 cm, and the dipole gaps vary from 14 cm to 9 cm. All 
magnet pole tip fields are under 2 T, except in the switchyard where they are as high as 
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TABLE X: Parameters for superconducting cavities. 





A 


B 


Frequency (MHz) 


201.25 


201.25 


Cells per cavity 


2 


2 


Aperture diameter (cm) 


46 


30 


Energy gain (MV) 


22.5 


25.5 


rf pulse length (ms) 


3 


3 


Input Power (kW) 


980 


1016 


Peak surface field (MV/m) 


23.1 


24.3 


Qo 


6 x 10 9 


6 x 10 9 



2.3 T in some cases (and the magnet apertures rise to 21 cm). 

Y-AA_/V\_/V\_AA^/, ,. 

FIG. 23: (Color)Layout of an RLA linac period. 

The 4.5 K-equivalent cryogenic load for the entire acceleration system is 27.9 kW. In 
Study I, where the final beam energy was chosen to be 50 GeV, a second RLA is needed. 
This second RLA is similar to the RLA just described, but considerably larger. 

G. Storage Ring 

After acceleration in the RLA, the muons are injected into the upward-going straight 
section of a racetrack-shaped storage ring with a circumference of 358 m. Parameters of 
the ring are summarized in Table IXll High-field superconducting arc magnets are used to 
minimize the arc length. Minimizing the arc length for a given length of straight maximizes 
the fraction of the circumference contained in the straight section, thereby maximizing the 
fraction of neutrinos (around 35% in our case) decaying toward the detector. 

Furthermore, the beta functions in the downward-going straight (which is pointed toward 
the detector) are made large to reduce the angular divergence of the beam. This ensures 
that the angular divergence of the beam is dominated by the calculable relativistic angular 
divergence of the decay neutrinos. The goal of this is not only to make the angular divergence 
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FIG. 24: (Color)Three-dimensional view of the storage ring magnets. 



of the neutrino beam as small as possible and therefore maximize the flux, but it also reduces 
the experimental uncertainty associated with an uncertainty in the flux that would come 
from an uncertainty in the angular divergence of the muon beam. 

All muons are allowed to decay; the maximum heat load from their decay electrons is 42 
kW (126 W/m). This load is too high to be dissipated in the superconducting coils. For 
Study II, a magnet design (see Fig. l2~4~)l has been chosen that allows the majority of these 
electrons to exit between separate upper and lower cryostats, and be dissipated in a dump at 
room temperature. To maintain the vertical cryostat separation in focusing elements, skew 
quadrupoles are employed in place of standard quadrupoles. The result is a skew FODO 
lattice, giving diagonal oscillations, as opposed to the horizontal and vertical oscillations of 
the usual upright FODO lattice. In order to maximize the average bending field, Nb 3 Sn 
pancake coils are employed. One coil of the bending magnet is extended and used as one 
half of the previous (or following) skew quadrupole to minimize unused space. Figure ESI 
shows a layout of the ring as it would be located at BNL. (The existing 110 m high BNL 
smokestack is shown for scale.) For site-specific reasons, the ring is kept above the local 
water table and is placed on a roughly 30 m high berm. This requirement places a premium 
on a compact storage ring. 

The beam is allowed to debunch in the storage ring. In one muon lifetime (0.42 ms), 
a bunch with the full energy spread (±2.2%) will have its total length increase by 0.51 /is 
(the storage ring is 1.19 /j,s long, and the bunch train starts out 0.33 /as long). If one 
wishes to avoid the increase in the bunch train length, one could install rf cavities, but the 
voltage required to avoid energy spread increase for the momentum compaction in this ring 
is prohibitive: a better solution would be a ring re-designed to have very low momentum 
compaction. 

For Study I, a conventional superconducting ring was utilized to store the 50 GeV muon 
beam. The heat load from muon decay products in this scenario is managed by having a liner 
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TABLE XI: Muon storage ring parameters. 



Energy (GeV) 


20 


Circumference (m) 


358.18 


Normalized transverse acceptance (mm rad) 


30 


Energy acceptance (%) 


2.2 


Momentum compaction 


0.028 


Arc 


Length (m) 


53.09 


No. cells per arc 


10 


Cell length (m) 


5.3 


Phase advance (deg) 


60 


Dipole length (m) 


1.89 


Dipole field (T) 


6.93 


Skew quadrupole length (m) 


0.76 


Skew quadrupole gradient (T/m) 


35 


Pmax (m) 


8.6 


Production Straight 


Length (m) 


126 


Anax (m) 


200 



inside the magnet bore to absorb the decay products. This approach is likewise available 
for BNL, provided some care is taken to keep the ring compact; acceptable lattice solutions 
have been found for this option as well. 

H. Overall Machine Summary 

Figure 123 shows the muons per incident proton in the "front end" (everything before the 
acceleration) of the Study II neutrino factory. Table IXHI gives the values at the ends of 
several sections and the losses in those sections. These significant losses are a necessary cost 
of making a low-emittance beam that can be accelerated and injected into a storage ring. 
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FIG. 25: Muons per incident proton in the Study II neutrino factory front end. 

TABLE XII: Muon survival and losses in several sections of the Study II neutrino factory. /j,/p is 
the number of muons per proton at the end of that section, and "Loss" is the loss in that section. 



Section 


m/p 


Loss 


Phase rotation 


0.60 




Buncher 


0.47 


22% 


Cooling 


0.22 


53% 


Accelerator aperture 


0.16 


26% 


Preaccelerator Linac 


0.15 


10% 


Recirculating Linac 


0.13 


10% 



An overall layout of the Neutrino Factory on the BNL site is shown in Fig. E3 Figure OKI 
shows the equivalent picture for a facility on the Fermilab site. In this latter case, the layout 
includes the additional RLA and longer storage ring needed to reach 50 GeV. Clearly the 
footprint of a Neutrino Factory is reasonably small, and such a machine would fit easily on 
the site of either BNL or Fermilab. 

I. Detector 

The Neutrino Factory, plus its long-baseline detector, will have a physics program that 
is a logical continuation of current and near-future neutrino oscillation experiments in the 
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FIG. 26: (Color)Top view and cross section through 20 GeV ring and berm. The existing 110 m 
tower, drawn to scale, gives a sense of the height of the ring on the BNL landscape. 

U.S., Japan, and Europe. Moreover, detector facilities located in experimental areas near 
the neutrino source will have access to integrated neutrino intensities 10 4 -10 5 times larger 
than previously available (10 20 neutrinos per year compared with 10 15 -10 16 ). 

The detector site taken for Study II is the Waste Isolation Pilot Plant (WIPP) in Carlsbad, 
New Mexico. The WIPP site is approximately 2900 km from BNL. Space is potentially 
available for a large underground physics facility at depths of 740-1100 m, and discussions 
are under way between DOE and the UNO project j^l on the possible development of such 
a facility. 

1. Far Detector 

Specifications for the long-baseline Neutrino Factory detector are rather typical for an 
accelerator-based neutrino experiment. However, the need to maintain a high neutrino rate 
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FIG. 27: Schematic of a 20 GeV Neutrino Factory at BNL. 



at these long distances requires detectors 3-10 times more massive than those in current 
neutrino experiments. Clearly, the rate of detected neutrinos depends on two factors — the 
source intensity and the detector size. In the final design of a Neutrino Factory, these two 
factors would be optimized together. 

Two options are considered for the WIPP site: a 50 kton steel-scintillator-proportional- 
drift-tube (PDT) detector or a water-Cherenkov detector. The PDT detector would resemble 
MINOS. Figure EE] shows a 50 kton detector with dimensions 8 m x 8 m x 150 m. A detector 
of this size would record up to 4 x 10 4 events per year. 

A large water-Cherenkov detector would be similar to SuperKamiokande, but with either 
a magnetized water volume or toroids separating smaller water tanks. The detector could be 
the UNO detector j^] , currently proposed to study both proton decay and cosmic neutrinos. 
UNO would be a 650 kton water-Cherenkov detector segmented into a minimum of three 
tanks (see Fig. I3T)|) . It would have an active fiducial mass of 440 kton and would record up 
to 3 x 10 5 events per year from the Neutrino Factory beam. 

Anot her p ossibility for a Neutrino Factory detector is a massive liquid-argon magnetized 



detector 
neutrinos 



that would also attempt to detect proton decay, as well as solar and supernova 



66 




FIG. 28: (Color) Schematic of a 50 GeV Neutrino Factory at Fermilab. 
2. Near Detector 

As noted, detector facilities located on-site at the Neutrino Factory would have access 
to unprecedented intensities of pure neutrino beams. This would enable standard neutrino 
physics studies, such as sin 2 ^, structure functions, neutrino cross sections, nuclear shad- 
owing and pQCD to be performed with much higher precision than previously obtainable. 
In addition to its primary physics program, the near detector can also provide a precise flux 
calibration for the far detector, though this may not be critical given the ability to monitor 
the storage ring beam intensity independently. 

A compact liquid-argon time projection chamber (TPC; similar to the ICARUS detec- 
tor [32}]), cylindrically shaped with a radius of 0.5 m and a length of 1 m, would have an 
active mass of 10 3 kg and a neutrino event rate O(10 Hz). The TPC could be combined 
with a downstream magnetic spectrometer for muon and hadron momentum measurements. 
At these neutrino intensities, it is even possible to have an experiment with a relatively thin 
Pb target (1 L ra d), followed by a standard fixed-target spectrometer containing tracking 
chambers, time-of-fiight, and calorimetry, with an event rate 0(1 Hz). 
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FIG. 29: (Color)A possible 50 kton steel-scintillator-PDT detector at WIPP. 




FIG. 30: (Color)Block schematic of the UNO detector, including initial design parameters. 
J. Staging Options 

It seems quite possible — perhaps even likely — that the Neutrino Factory would be built 
in stages, both for programmatic and for cost reasons. Here we outline a possible staging 
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concept that provides good physics opportunities at each stage. The staging scenario we 
consider is not unique, nor is it necessarily optimized. Depending on the results of our 
technical studies and the results of continued searches for the Higgs boson, it is hoped that 
the Neutrino Factory is really the penultimate stage, to be followed later by a Muon Collider 
(Higgs Factory). We assume this possibility in the staging discussion that follows. Because 
the physics program would be different at different stages, it is impractical at this time to 
consider detector details. 



1. Stage 1 

In the first stage, we envision a Proton Driver and a Target Facility to create superbeams. 
The Driver could begin with a 1 MW beam level (Stage 1) or could be designed from the 
outset to reach 4 MW (Stage la). (Since the cost differential between 1 and 4 MW is not 
expected to be large, we do not consider any intermediate options here.) It is assumed, as 
was the case for both Study I and Study II, that the Target Facility is built from the outset 
to accommodate a 4 MW beam. Based on the Study II results, a 1 MW beam would provide 
about 1.2 x 10 14 /i/s (1.2 x 10 21 /i/year) and a 4 MW beam about 5 x 10 14 /i/s (5 x 10 21 
/i/year) into a solenoid channel. 

In addition t o the neu trino program, this stage will also benefit it, K, and p programs, 



as discussed in 
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2. Stage 2 

In Stage 2, we envision a muon beam that has been phase rotated (to give a reasonably 
low momentum spread) and transversely cooled. In the nomenclature of Study II, this stage 
takes us to the end of the cooling channel. Thus, we have access to a muon beam with a 
central momentum of about 200 MeV/c, a transverse (normalized) emittance of 2.7 mm rad 
and an rms energy spread of about 4.5%. The intensity of the beam would be about 4 x 10 13 
/i/s (4 x 10 20 /i/year) at 1 MW, or 1.7 x 10 14 /i/s (1.7 x 10 21 /i/year) at 4 MW. If more 
intensity were needed, and if less cooling could be tolerated, the length of the cooling channel 
could be reduced. As an example, stopping at the end of Lattice 1 instead of the end of 
Lattice 2 in the Study II cooling channel would result in an increase of transverse emittance 
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by roughly a factor of two. This is an appropriate stage to mount an experiment to search 
for a non-zero muon electric dipole moment. 

3. Stage 3 

In Stage 3, we envision using the Pre-acceleration Linac to raise the beam energy to 
roughly 2.5 GeV. At this juncture, it may be appropriate to consider a small storage ring, 
comparable to the g— 2 ring at BNL, to be used for the next round of muon g — 2 experiments. 

4- Stage 4 

At Stage 4, we envision having a complete Neutrino Factory operating with a 20 GeV 
storage ring. This stage includes the RLA and the storage ring. If it were necessary to 
provide a 50 GeV muon beam as Stage 4a, an additional RLA and a larger storage ring 
would be needed. 

5. Stage 5 

In Stage 5, we could envision an entry-level Muon Collider operating as a Higgs Factory. 
Because the initial muon beam must be prepared as a single bunch of each charge, an addi- 
tional ring for the proton driver to coalesce proton bunches into a single pulse is anticipated. 
The cooling will have to be significantly augmented. First, a much lower transverse emit- 
tance is needed, and second, it will be necessary to provide longitudinal cooling (emittance 
exchange) to maintain a reasonable transmission of the muons. The additional cooling will 
permit going to smaller solenoids and higher frequency rf systems (402.5 or perhaps 805 
MHz), which should provide more cost-effective cooling. Next, we will need considerably 
more acceleration, though with smaller energy acceptance and aperture requirements than 
at present. Lastly, we will need a very low (3* lattice for the collider ring, along with miti- 
gation of the potentially copious background levels near the interaction point. In this case 
the detector is, in effect, part of the collider ring, and its design must be an integral part of 
the overall ring design. 
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IV. MUON COLLIDERS 



The primary advantage of using muons in a lepton collider arises from the fact that the 
muon is ~ 200 times heavier than the electron. It is thus possible to accelerate muons using 
circular accelerators that are compact and fit on existing accelerator sites. See Figure 
for a comparison of relative sizes of muon colliders ranging from 500 GeV to 3 TeV center 
of mass energies with respect to the LHC, SSC, and NLC. Once the problem of cooling 
a muon beam to sufficiently small emittances is solved and the beams can be accelerated, 
higher energies are much more easily obtained in a muon collider than in the linear electron- 
positron collider. Because the muon is unstable, it is necessary to cool and accelerate the 
beam before a substantial number have decayed. The number of turns in a muon lifetime is 
independent of the muon momentum for a given magnetic field, since both the circumference 
and the muon lifetime in the laboratory frame scale with muon momentum. With typical 
bending magnetic fields(~ 5 Tesla) available with today's technology, the muons last « 1000 
turns before half of them have decayed in the collider ring. 

Muon decay also gives rise to large numbers of electrons that can affect the cryogenics of 
the magnets in the machine as well as pose serious background problems for detectors in the 
collision region. The 1999 Status Report Q] contains an excellent summary of the problems 
(and possible solutions) one faces on the way to a muon collider. 

Figure E21 shows a schematic of such a muon collider, along with a depiction of the possible 
physics that can be addressed with each stage of the facility. Some of the ideas needed to 
obtain longitudinal cooling necessary for the Muon Collider are discussed in section II V Bl and 
some of the parameters of the accelerator system for higher energy colliders are discussed in 
section IIVCI below. 



A. Higgs Factory Requirements 

The emittance of the muon beam needs to be reduced by a factor of rs 10 6 from pro- 
duction JqJ to the point of collision for there to be significant luminosity for experiments. 
Table IXIIII li sts the transverse and longitudinal emittances at the end of the decay channel, 
Study-II [29! cooling channel and those needed for a 0.1 TeV Center of Mass Energy Muon 
Collider, also known as a Higgs Factory. It can be seen that one needs to cool by a factor of 
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(5-17 TeV) 



FIG. 31: (Color) Comparative sizes of various proposed high energy colliders compared with the 
FNAL and BNL sites. The energies in parentheses give for lepton colliders their CoM energies 
and for hadron colliders the approximate range of CoM energies attainable for hard parton-parton 
collisions. 

« 20 in the transverse dimension and ~ 3 in the longitudinal dimension from the Study-II 
emittances to achieve the emittances necessary for a Higgs factory. This can be achieved by 

TABLE XIII: Transverse and longitudinal emittances at the end of the decay channel, Study-II 



5 channel and the Hig 


gs factory cooling channel. 




Emittance at end of 


Transverse emittance (ir mm) 


Longitudinal emittance (ir mm) 


Decay Channel 


17 


150 


Study-II cooler 


2.6 


30 


Higgs Factory cooler 


0.14 


9 
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FIG. 32: (Color)Schematic of a muon collider. 

ionization cooling similar to the scheme described in the section IIHI The transverse emit- 
tance is reduced during ionization cooling, since only the longitudinal energy loss is replaced 
by rf acceleration. However, due to straggling, the longitudinal energy spread of the beam 
increases, even if the average longitudinal energy of the beam is kept constant. The longitu- 
dinal emittance thus grows in a linear cooling channel. In order to cool longitudinally, one 
needs to create dispersion in the system and have wedge absorbers at the point of maximum 
dispersion so that the faster particles go through the thicker parts of the wedge. This results 
in a reduction in longitudinal emittance accompanied by an increase in transverse emittance 
and is thus called emittance exchange. 

The status r epo rt Q outhnes the detaHs of the aeration and eo,hde r n„ g for the 
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0.1 TeV Higgs factory, shown schematically in Fig. ESI Table IXIVl gives a summary of the 
parameters of various muon colliders including three different modes of running the Higgs 
Collider that have varying bea m momen tum spreads. Additional information about the 



Muon Collider can be found in 
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TABLE XIV: Baseline parameters for high- and low-energy muon colliders. Higgs/year assumes a 
cross section a = 5 x 10 4 fb; a Higgs width T = 2.7 MeV; 1 year = 10 s. 
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FIG. 33: (Color)Plan of a 0.1-TeV-CoM muon collider, also known as a Higgs Factory. 
B. Longitudinal Cooling 

At the time of writing of the status report [^] there was no satisfactory solution for the 
emittance exchange problem and this remained a major stumbling block towards realizing a 
muon collider. However, ring coolers have been found to hold significant promise in cooling 
in 6-dimensional phase space. Another advantage of ring coolers is that one can circulate 
the muons many turns, thereb y re using the cooling channel el emen ts. Several meetings on 
emittance exchange were held jl4l| and a successful workshop |l42j | was held in 2001, where 
we explored in some depth several kinds of ring coolers. These options differ primarily 
in the type of focusing used to contain the beam. We describe the current status of our 
understanding of three types of ring coolers here. 



1. Solenoidal Ring Coolers 



The basic design of the solenoidal ring cooler |143| is presented in Figure 0U Eight 
focusing dipole magnets with an index n = — ~ are used for bending and focusing of the 
beam. Each of these dipoles bends the beam through 45 degrees with a central orbit bending 
radius of 52 cm. We have done calculations to show that such dipoles are buildable. FigurelHol 
shows a configuration of such a dipole and the resulting magnetic field components calculated 
using a 3D field calculation program. There are 4 long solenoids containing rf cavities and 
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liquid hydrogen absorbers for transverse cooling. A magnetic field of 2.06 T at the end 
regions of the solenoids provides the same transverse focusing as the bending magnets. The 
magnetic field adiabatically increases to 5.15 T towards the center of the solenoid in order 
to produce a small (3 function (25-30 cm) at the absorbers. The short solenoids are designed 
to create an appropriate dispersion function that is zero at the long solenoids, which house 
the 200 MHz rf cavities. Their field is ±2.06 T at the edges and ±2.75 T centrally. A 
symmetric field flip is required in the short solenoids to prevent the build up of canonical 
angular momentum. This field flip causes the dispersion in the long solenoids housing the 
rf cavities to be zero while permitting a non-zero dispersion at the lithium hydride wedge 
absorbers at the centers of the short solenoids which then produce longitudinal cooling via 
emittance exchange. 
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FIG. 34: (Color)Layout and parameters of the solenoid based ring cooler 



Evolution of the beam emittance and transmission is shown in Figure EHl a s a function 
of the number of turns in the ring. In 15 turns, the transverse emittance decreases from 1.2 
cm to 0.21 cm yielding a cooling factor of 5.7, the longitudinal emittance decreases from 1.5 
cm to 0.63 cm (cooling factor 2.4), and the 6-D emittance decreases from 2.2 cm 3 to 0.028 
cm 3 , with an overall cooling factor « 79. The transmission is 0.71 without decay and 0.48 
with decay. We define a merit factor for cooling that is the total transmission including 
decay times the 6-D cooling factor. The merit factor for this ring is then 38. This implies 
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FIG. 35: (Color)Evolution of the beam emittance/transmission at the ring cooler. 
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FIG. 36: (Color)Figure shows a computer model of a 52cm radius dipole with index n = — | and 
the calculated field components radially 10cm off axis. 



that transverse emittance at the ring cooler is about the same as at a linear SFOFO cooling 
channel employed in Study-II Q], whereas the longitudinal emittance is noticeably less. 

This cooler provides mainly transverse cooling and can be used as a part of Neutrino 
Factory or a muon collider. A cooler specially designed for strong longitudinal cooling 
( "bunch compressor" ) can also be created using a similar scheme. Such a compressor would 
be a part of a muon collider to shorten muon bunches from 6-8 m (minimal length after 
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3: 



n — jj decay and phase rotation, see Ref. |9|) to 0.6-0.8 m acceptable for further cooling by 
a 200 MHz channel. 

Two options for the bunch compressor are considered in Ref. 1441 ]. The first one is a two- 
step cooler where each step is very similar to the ring cooler shown in Figure The main 
difference is that the primary goal in the first cooler is the longitudinal bunching of the beam. 
This leads to a uniform magnetic field in the long solenoids and lower frequency/voltage of 
the accelerating rf system (15.6 MHz/4 MeV/m at the first stage vs. 62.5 MHz/8 MeV/m 
at the second one). Another option is a 15 MHz octagonal cooler composed of the same cells 
as in Figure EU but with half the bending magnet angle. Decrease of longitudinal emittance 
from 43 cm to 2.5-3 cm, as required for muon collider, is obtained in both cases. 

We are proceeding with a realistic sim ulation of this system using Geant and ICOOL 
that employs realistic magnetic fields jl4sl | produced by field calculation programs. 

After the two stage cooler, we still need a factor of ~ 30 in transverse cooling, but we 
are within a factor of 4 in longitudinal cooling relative to the Higgs factory goals. Lithium 
lens cooling, which with its strong focusing will cool transversely further while degrading 
longitudinally due to straggling, is a possibility and is being investigated. 



2. RFOFO ring coolers 

The cooling lattice for the Neutrino Factory (see section IHIjl employs a configuration 
of fields known as an SFOFO lattice (super-FOFO) where the axial magnetic field profile 
changes polarity in alternate cells of the lattice. For the ring cooler design under considera- 
tion here, we employ an RFOFO lattice (regular-FOFO) where the axial field profile changes 
polarity in the middle of a cell. As a result all cells in an RFOFO lattice are identical. 

The ring cooler design employs a single cell for both transverse cooling and emittance 
exchange. It uses solenoids for focusing, giving large angular and momentum acceptances. 
The cell includes dispersion, acceleration, and energy loss in a single thick hydrogen wedge. 
Figure shows the layout of the cooling ring drawn to scale. The RFOFO lattice was chosen 
because, unlike in the SFOFO case used in Study-II, all cells are strictly identical, and the 
presence of an integer betatron resonance within the momentum acceptance is eliminated. 

The basic 33 m circumference ring is made up of 12 identical 2.75-m long cells. In the 
figure, this symmetry is broken for injection and extraction, but the magnetic fields in this 
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FIG. 37: (Color)Layout of an RFOFO cooling ring. 

insertion are nearly identical to those in the rest of the ring. Figure |58] shows a detailed 
view of three cells of the lattice, in plan (a) and side (b) views. 

The longitudinal field on-axis has an approximately sinusoidal dependence on position. 
The beam axis is displaced laterally with respect to the coil centers (as shown in figure 
158k ) to minimize horizontal fields that cause vertical beam deviations. The lattice transmits 
particles in the momentum band from 150 to 250 MeV/c. The minimum value of the beta 
function at the central momentum is 40 cm. 

Bending round the ring, and the required dispersion, are provided by applying an ap- 
proximately 0.125 T vertical bending field generated by alternately tilting the solenoids (as 
shown in figure 188b). There is no attempt to set the field index n (where B oc r n ) to the 
0.5 value; so the focusing is not identical in x and y. 

It is found that the acceptance is reduced as the bending field is increased. We thus use 
a wedge with maximum possible angle (giving zero thickness on one side), and the lowest 
bending field consistent with adequate emittance exchange. The dispersion at the absorber 
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FIG. 38: (Color)Three cells of the RFOFO lattice; a) plan, b) side. 

is approximately 8 cm in a direction 30 degrees from the y axis, The dispersion at the center 
of the rf is of the opposite sign, and also mostly in the y direction. 

The liquid-hydrogen wedge has a central thickness of 28.6 cm and a total wedge angle of 
100 degrees and is rotated 30 degrees from the vertical. The rf cavities are at a frequency 



of 201.25 MHz and 
The ICOOL 



rave an accelerating gradient of 12 MV/m. 
145l ] simulation (with results shown in shown Figure l3*^j) used fields calcu- 
lated from the specified coils, and thus neglects no end field effects. But in this simulation, 
no absorber, or rf, windows, are included; nor did it include the injection/extraction in- 
sertion. The rf was represented by the fields in perfect pillbox cavities. The input tracks 
were taken from a Study-II 2f| simulation, using distributions from just upstream of the 
transverse cooling system. The use of Study-II simulated distributions is intended to allow 
a more realistic estimate of the ring's performance. No attempt was made to match this 
beam to the ring dispersion or slight differences in the transverse beta functions. 

Figure EH shows the simulated transmission, transverse emittance, longitudinal emittance, 
6-dimensional emittance, and a merit factor M vs. length in the ring. M is given by: 

M = 66 V 7~ i\ x Transmission 
e 6 (fznal) 

Initially, the x emittance falls more rapidly than the y, because it is the y emittance 
that is exchanged with the longitudinal. But the Larmor rotations soon mix the x and y 
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FIG. 39: (Color)Transmission, normalized transverse emittance, normalized longitudinal emit- 
tance, normalized 6-dimensional emittance, and the merit factor, as a function of distance. 

emittances bringing them to a common value. 

After a distance of 400 m (« 12 turns), the 6-dimensional emittance has fallen by a factor 
of 238, with a transmission of 48 % (66% without decay). The merit factor is 136. The same 
factor for the Study-II cooling lattice, also without windows, is 13. Studies with realistic 
windows and the injection/extraction insertion added, show lower merit factors, but always 
far better than the Study-II example. 

The design of the injection/ejection channels and kickers will be challenging, and this 
ring could not be used, as is, to replace the Study-II cooling channel because the bunch 
train in that case is too long to fit in the ring. Both problems would be removed in a helical 
cooling channel. The merit factor for such a channel could be even better than that of the 
ring because it would be possible to "taper" the optics, as a function of distance down the 
channel, and thus lower the final equilibrium emittance. 
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3. Quadrupole Ring Coolers 



Alternative ring designs have been explored which are based on storage rings w hich c onsist 



of conventional quadrupole and dipole magnetic elements instead of sol enoids |149j . The 
strategy has been to utilize the SYNCH storage ring design code 152j ] to develop linear 
lattice s oluti ons and then transfer the lattice parameters into the ray-tracing tracking code 



ICOOL 



145| in which absorbers and energy recovery with rf cavities can be included for full 



simulation. An example of such a ring is shown in Figure EHl where the elements of a half 
cell for a 22.5° bending cell are depicted schematically. The correspondence between the 
beam envelope beta and dispersion functions resulting from a simulation with the ICOOL 
tracking code and the SYNCH calculated values are shown in Figure |^ The full sixteen 
cell ring is shown in Figure | 
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FIG. 40: (Color) Schematic diagram of half of a 22.5° bending cell. A wedge absorber is located 
in the middle of the cell. 



In general, we find that the performance of the rings, as measured in emittance reduction 
along with particle transmission and decay losses, improves when more compact lattice 
designs are considered. In Figure 021 the variation of the normalized emittances as a function 
of ring turns is shown for an eight cell ring. A reduction of normalized emittance is observed 
for all three dimensions. This particular ring has a total circumference of 30.9 m. Each 
half cell contains one 22.5° combined function dipole proceeded and followed by a single 
horizontally focusing quadrupole. The average muon beam momentum is 250 MeV/c and 
liquid hydrogen absorbers with wedge opening angles of 40° are used. For each cell, the 



82 



P x , p y , lO* D <m) 22.5 deg. Bent Lattice 





1 1 II 


1 1 VI 


1 III 


1 


J 


: rf 


Q B 


Q Q 
Wedge 


B O 


RF 




absorber 
Marked pointS(ICOOL) : 






SiOl 


a curvc5(»YNun) 








- 0* 




























• • 




























10* 


D 











s <m) 



FIG. 41: (Color)The and D (dispersion) in a 22.5° bending cell. SYNCH calculations (solid 

curves) and beam parameters from an ICOOL simulation (marked points) are compared. 




FIG. 42: (Color) Top view of a sixteen cell muon cooling ring. 

central beam orbit traverses 24 cm of absorber. The energy loss in the wedge absorbers is 
compensated with 201 MHz rf cavities with peak on-axis gradients of 16 MV/m. 

The muon transmission efficiency and total merit factor (muon survival rate times the 
ratio of initial to final 6-dimensional emittance) is shown in Figure EU as a function of ring 
turns. The merit factor reaches 18, while the muon transmission efficiency including decay 
losses after eighteen full turns is 40%. 

Rings in which the focusing function is handled exclusively by the dipole elements have 
also been explored. In this case, the natural focusing power of the dipole is utilized for hor- 
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FIG. 43: (Color)The evolution of x, y, z normalized emittances in 32 full turns. 
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FIG. 44: (Color) The transmission and the figure of merit factor as a function of the arc length in 
32 full turns. 

izontal focusing while the entrance and exit dipole edge angles are adjusted to provide the 
required vertical focusing. Several examples of such lattices have been examined. Because 
these lattices are more compact than lattices which also include quadrupoles, the perfor- 
mance of these rings tend to be better. As an example, performances with merit factors 
on the order of 100 have been observed with rings based on four cell lattice with eight 45° 
dipoles. For one case, the entrance and exit angles dipole faces are 7.4° and 21.8° relative 
to the normal of the beam trajectory. The ring for this example has a circumference of only 
9.8 m and the design of injection/ejection cells may prove to be very challenging. 
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4- Injection into Ring Coolers 



The most serious technical problem facing the ri ng co oler approach is the injection system 



which may require a very powerful kicker magnet |146| . The energy stored in the injection 
kicker goes as the square of the emittance of the beam and inversely as the circumference of 
the ring. A promising injection scheme that does not use kicker magnets, but instead uses 



absorbers to degrade the beam energy and rf phase manipulations has been proposed 
and is being studied. 



147| 



C. Higher Energy Muon colliders 

Once the required cooling has been achieved to make the first muon collider feasible, 
acceleration to higher energies becomes possible. Colliders with 4 TeV center of mass energy 
have been studied P] and Table IX VI lists the parameters for such a collider. The muons 
are accelerated initially by a linear accelerator followed by a series of recirculating linear 
accelerators (RLA's) followed by rapid cycling synchrotrons (RCS's). The radiation from 
the neutr inos from the muon decay begins to become a problem at CoM energies of 3 
TeV Q 

There have been preliminary attempts to study colliders of even higher energy, starting 
at 10 TeV all the way up to 100 TeV in the center of mass and we include the references to 



these studies 



1531 ] for the sake of completeness. 



D. Muon Collider Detectors 

Figure 021 shows a strawman muon collider detector for a Higgs factory simulated in 
Geant. The background from muon decay sources has been extensively studied JjJ. At the 
Higgs factory, the main sources of background are from photons generated by the showering 
of muon decay electrons. At the higher energy colliders, Bethe-Heitler muons produced in 
electron showers become a problem. Work was done to optimize the shielding by using spe- 
cially shaped tungsten cones [l| that reduce the backgrounds resulting from electomagnetic 
showers from entering the detector. The occupancy levels due to background photons and 
neutrons in detectors were shown to be similar to those predicted for the LHC experiments. 
It still needs to be established whether pattern recognition is possible in the presence of these 
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TABLE XV: Parameters of Acceleration for a 4 TeV Muon Collider. 





Linac 


RLA1 


RLA2 


RCS1 


RCS2 


E (GeV) 


0.1-> 1.5 


1.5 10 


10 70 


70 250 


250 2000 


frf (MHz) 


30 100 


200 


400 


800 


1300 


^turns 


1 


9 


11 


33 


45 


V r/ (GV/turn) 


1.5 


1.0 


6 


6.5 


42 


C 4urn (km) 


0.3 


0.16 


1.1 


2.0 


11.5 


Beam transit time (ms) 


0.0013 


0.005 


0.04 


0.22 


1.73 


Vz,beam (™) 


50 -> 8 


4 -> 1.7 


1.7 0.5 


0.5 0.25 


0.25 -> 0.12 




0.005 -> 0.033 


0.067 0.16 


0.16 0.58 


0.58 1.14 


1.14 -> 2.3 


Loss (%) 


5 


7 


6 


7 


10 



backgrounds, especially the Bethe Heitler muons, which are a unique source of background 
to muon collider detectors.. 




FIG. 45: (Color)Cut view of a strawman detector in Geant for the Higgs factory with a Higgs— > bb 
event superimposed. No backgrounds shown. The tungsten cones on either side of the interaction 
region mask out a 20 deg area. 
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V. R&D STATUS 



A successful construction of a muon storage ring to provide a copious source of neutrinos 
requires several novel approaches to be developed and demonstrated; a high-luminosity Muon 
Collider requires an even greater extension of the present state of accelerator design. Thus, 
reaching the full facility performance in either case requires an extensive R&D program. 

The required R&D program has been identified for each of the major neutrino factory 
systems. In particular, some critical components must be prototyped and their performance 
verified. For example, the cooling channel assumes a normal conducting rf (NCRF) cavity 
gradient of 17 MV/m at 201.25 MHz, and the acceleration section demands similar per- 
formance from superconducting rf (SCRF) cavities at this frequency. In both cases, the 
requirements are beyond the performance reached to date for cavities in this frequency 
range. The ability of the target to withstand a proton beam power of up to 4 MW must 
be confirmed, and, if it remains the technology of choice, the ability of an induction linac 
unit to coexist with its internal SC solenoid must be verified. Finally, an ionization cooling 
experiment should be undertaken to validate the implementation and performance of the 
cooling channel, and to confirm that our simulations of the cooling process are accurate. 

Below we give an overview of the MC R&D program goals and list the specific questions 
we expect address. We also summarize briefly the R&D accomplishments to date 

A. R&D Program Overview 

A Neutrino Factory comprises the following major systems: Proton Driver, Target and 
(Pion) Capture Section, (Pion-to-Muon) Decay and Phase Rotation Section, Bunching and 
Matching Section, Cooling Section, Acceleration Section, and Storage Ring. These same 
categories apply to a Muon Collider, with the exception that the Storage Ring is replaced 
by a Collider Ring having a low-beta interaction point and a local detector. Parameters 
and requirements for the various systems are generally more severe in the case of the Muon 
Collider, so a Neutrino Factory can properly be viewed as a scientifically productive first 
step toward the eventual goal of a collider. 

The R&D program we envision is designed to answer the key questions needed to embark 
upon a Zeroth-order Design Report (ZDR). The ZDR will examine the complete set of 
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systems needed for a Neutrino Factory, and show how the parts can be integrated into a 
coherent whole. Although a fully engineered design with a detailed cost estimate is beyond 
the scope of a ZDR, enough detailed work must been accomplished to ensure that the critical 
items are technically feasible and that the proposed facility could be successfully constructed 
and operated at its design specifications. By the end of the full R&D program, it is expected 
that a formal Conceptual Design Report (CDR) for a Neutrino Factory could be written. 
The CDR would document a complete and fully engineered design for the facility, including 
a detailed bottom-up cost estimate for all components. This document would form the basis 
for a full technical, cost, and schedule review of the construction proposal. Construction 
could then commence after obtaining approval. 

The R&D issues for each of the major systems must be addressed by a mix of theoretical 
calculation, simulation modeling, and experimental studies, as appropriate. A list of the key 
physics and technology issues for each major system is given:- 

Proton Driver 

• Production of intense, short proton bunches, e.g., with space-charge compensation 
and/or high-gradient, low frequency rf systems 

Target and Capture Section 

• Optimization of target material (low-Z or high-Z) and form (solid, moving band, 
liquid-metal jet) 

• Design and performance of a high-field solenoid (^20 T) in a very high radiation 
environment 

Decay and Phase Rotation Section 

• Development of high-gradient induction linac modules having an internal supercon- 
ducting solenoid channel 

• Examination of alternative approaches, e.g., based upon combined rf phase rotation 
and bunching systems or fixed-field, alternating gradient (FFAG) rings 

Bunching and Matching Section 
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• Design of efficient and cost-effective bunching system 
Cooling Section 

• Development and testing of high-gradient normal conducting rf (NCRF) cavities at a 
frequency near 200 MHz 

• Development and testing of efficient high-power rf sources at a frequency near 200 
MHz 

• Development and testing of liquid hydrogen absorbers for muon cooling 

• Development and testing of an alternative gaseous-absorber cooling-channel design 
incorporating pressurized, high-gradient rf cavities. 

• Development and testing of candidate diagnostics to measure emittance and optimize 
cooling channel performance 

• Design of beamline and test setup (e.g., detectors) needed for demonstration of trans- 
verse emittance cooling 

• Development of full six-dimensional analytical theory to guide the design of the cooling 
section 

Acceleration Section 

• Optimization of acceleration techniques to increase the energy of a muon beam (with 
a large momentum spread) from a few GeV to a few tens of GeV (e.g., recirculating 
linacs, rapid cycling synchrotrons, FFAG rings) for a Neutrino Factory, or even higher 
energy for a Muon Collider 

• Development of high-gradient superconducting rf (SCRF) cavities at frequencies near 
200 MHz, along with efficient power sources (about 10 MW peak) to drive them 

• Design and testing of components (rf cavities, magnets, diagnostics) that will operate 
in the muon-decay radiation environment 

Storage Ring 
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• Design of large- aperture, well-shielded superconducting magnets that will operate in 
the muon-decay radiation environment 

Collider 

• Cooling of 6D emittance (x, p x , y, p y , t, E) by up to a factor of 10 5 — 10 6 

• Design of a collider ring with very low (3* (a few mm) at the interaction point having 
sufficient dynamic aperture to maintain luminosity for about 500 turns 

• Study of muon beam dynamics at large longitudinal space-charge parameter and at 
high beam-beam tune shift 

Detector 

• A study of cost trade-offs between increasing the detector mass compared with in- 
creasing the beam intensity. 

• Simulation studies to define acceptable approaches for a Muon Collider detector op- 
erating in a high-background environment 

Most of these issues are being actively pursued as part of the ongoing MC R&D program. 
In a few areas, notably the proton driver and detector, the MC does not currently engage in 
R&D activities, though independent efforts are under way. Longer-term activities, related 
primarily to the Muon Collider, are supported and encouraged as resources permit. 

B. Recent R&D Accomplishments 

1. Targetry 

A primary effort of the Targetry experiment E951 at BNL has been to carry out initial 
beam tests of both a solid carbon target and a mercury target. Both of these tests have 
been made with a beam intensity of about 4 x 10 12 ppp, with encouraging results. 

In the case of the solid carbon target, it was found that a carbon-carbon composite having 
nearly zero coefficient of thermal expansion is largely immune to beam-induced pressure 
waves. Sublimation losses are a concern. However, a carbon target in a helium atmosphere 
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is expected to have negligible sublimation loss. If radiation damage is the limiting effect for 
a carbon target, the predicted lifetime would be about 12 weeks when bombarded with a 
1 MW proton beam. Carbon targets would therefore seem viable for beam powers up to 
1 MW, or perhaps a little higher. 

For a mercury jet target, tests with about 2 x 10 12 ppp showed that the jet is not dispersed 
until long after the beam pulse has passed through the target. Measurements of the velocity 
of droplets emanating from the jet as it is hit with the proton beam pulse from the AGS (~10 
m/s for 25 J/g energy deposition) compare favorably with simulation estimates. High-speed 
photographs indicate that the beam disruption at the present intensity does not propagate 
back upstream toward the jet nozzle. If this remains true at the higher intensity of 1.6 x 10 13 
ppp, it will ease mechanical design issues for the nozzle. The tests so far have used a jet with 
an initial mercury velocity of 2 m/sec. At a Neutrino Factory a 20 m/sec jet is envisioned. 
A prototype is under development. 



2. MUCOOL 



The primary effort in the muon ionization cooling (MUCOOL) R&D has been to complete 
the Lab G high power rf test area at Fermilab and begin high-power tests of 805-MHz rf 
cavities. A test solenoid for the facility, capable of operating either in solenoid mode (its 
two independent coils powered in the same polarity) or gradient mode (with the two coils 
opposed), was commissioned up to its design field of 5 T. 

An 805 MHz open-cell cavity has been tested in Lab G to look at gradient limitations, 
magnetic field effects and compatibility of the rf cavities with other systems. We have mea- 
sured the dark currents over a range covering 14 orders of magnitude, and accumulated data 
on the momentum spectrum, angular di strib ution, pulse shape, dependence on condition- 
ing and dependence on magnetic fields jl56 |. The dark currents seem to be described by 
the Fowler Nordheim field emission process, which results from very small emitter sources 
(sub-micron sizes) at very high local electric fields (5-8 GV/m). This implies that the 
emitter fields are enhanced by large factors, (3fn =~ 500, over the accelerating field. (At 
these electric fields the electrostatic stress becomes comparable to the strength of hardened 
copper.) We have shown how both normal conditioning and nitrogen processing can reduce 
dark currents. Our data from the 805 MHz cavity has been compared with other data from 



91 



NLC cavities, superconducting TESLA cavities and 200 MHz proton linacs, showing that 
all cavities seem to be affected by the same processes. 

We have also looked at damage produced on irises and windows, primarily when the 
system is run with the solenoid magnet on. A number of effects are seen: copper splatters 
on the inside of the thin Ti window, burn marks on the outside of the window due to electron 
beamlets, and some craters, evidently produced by breakdown on the irises. The electron 
beamlets burned through the windows twice. We have measured the parameters of the 
beamlets produced from individual emitters when the magnetic field is on, and we have seen 
ring beams, presumably produced by ExB drifts during the period when the electrons are 
being accelerated. The radius of the beamlets is found to be proportional to E/B 2 . 

We are proceeding with an experimental program designed to minimize the dark currents 
using surface treatment of the copper cavity. 

A second cavity, a single-cell pillbox having foils to close the beam iris, has been tuned 
to final frequency, shipped to Fermilab, and testing has begun at Lab G. This cavity will 
permit an assessment of the behavior of the foils under rf heating and give indications about 
multipactor effects. It will also be used to study the dark current effects discussed above. 
An advantage of the pillbox cavity is that its windows can be replaced with ones made from 
(or coated with) various materials and cleaned or polished by various techniques. 

Development of a prototype liquid hydrogen absorber is in progress. Several large diame- 
ter, thin (125 fj,m) aluminum windows have been successfully fabricated by machining from 
solid disks. These have been pressure tested with water and found to break at a pressure 
consistent with design calculations. A new area, the MUCOOL Test Area (MTA), is being 
developed at FNAL for testing the absorbers. The MTA, located at the end of the proton 
linac, will be designed to eventually permit beam tests of components and detectors with 
400 MeV protons. It will also have access to 201-MHz high-power rf amplifiers for testing 
of future full-sized 201-MHz cavities. 

Initial plans for a cooling demonstration are well under way. This topic is covered sepa- 
rately in more detail in Section IVfl 

A parallel cooling channel development effort base d on the use of gaseous hydrogen or 



154j | has received a DOE Small Business 



helium energy-absorber has begun. Muons Inc. 
Technology Transfer grant with the Illinois Institute of Technology to develop cold, pres- 
surized high-gradient rf cavities for use in muon ionization cooling. These cavities will be 
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filled with dense gas, which suppresses high voltage breakdown by virtue of the Paschen 
effect and also serves as the energy absorber. A program of development for this alternative 
approach to ionization cooling is foreseen that starts with Lab G tests, evolves to an MTA 
measurement program, and leads to the construction of a cooling channel section suitable 
for tests in MICE. 

3. Feasibility Study-II 

The MC has participated in a second Feasibility Study for a Neutrino Factory, co- 
sponsored by BNL. The results of the study were encouraging (see Section III), indicating 
that a neutrino intensity of 1 x 10 20 per Snowmass year per MW can be sent to a detec- 
tor located 3000 km from the muon storage ring. It was also clearly demonstrated that a 
Neutrino Factory could be sited at either FNAL or BNL. Hardware R&D needed for such a 
facility was identified, and is a major part of the program outlined here. 

4- Beam Simulations and Theory 

I n ad dition to work on Study-II, our present effort has focused on longitudinal dynam- 
ics jl55l |. We are developing theoretical tools for understanding the longitudinal aspects of 
cooling, with the goal of developing approaches to 6D cooling, i.e., "emittance exchange." 
This is a crucial aspect for the eventual development of a Muon Collider, and would also 
benefit a Neutrino Factory. Improved designs for the phase rotation, bunching, and acceler- 
ation systems are also being explored, with an emphasis on preparing the way for a future 
design study in which the performance obtained by the study II design is maintained, but 
with a reduction in cost. 

5. Other Component Development 

At present, the main effort in this area is aimed at development of a high-gradient 201- 
MHz SCRF cavity. A test area of suitable dimensions has been constructed at Cornell. 
In addition, a prototype cavity has been fabricated for the Cornell group by our CERN 
colleagues. Mechanical engineering studies of microphonics and Lorentz detuning issues are 
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being carried out. 



6. Collider R&D 

Studies of possible hardware conf igura tions to perform emittance exchange, such as the 
compact ring proposed by Balbekov |l43j . are now getting under way. A ring cooler has the 
potential to cool in 6D phase space, provided the beam can be injected into and extracted 
from it. A series of workshops have been held on the topics of emittance exchange and ring 
coolers that have helped further our understanding of both. 

VI. INTERNATIONAL MUON IONIZATION COOLING EXPERIMENT 
A. Motivation 

Ionization cooling of minimum-ionizing muons is an important ingredient in the perfor- 
mance of a Neutrino Factory. However, it has not been demonstrated experimentally. We 
seek to carry out an experimental demonstration of cooling in a muon beam. Towards this 
goal, we have developed (in collaboration with a number of physicists from Europe and 
Japan interested in neutrino factories) a conceptual design for an International Muon Ion- 
ization Cooling Experiment (MICE). A proposal for MICE has recently been submitted to 



24] 



the Rutherford Appleton Laboratory in England 

The aim of the proposed cooling experimental demonstration is 

• to show that we can design, engineer and build a section of cooling channel capable of 
giving the desired performance for a neutrino factory; 

• to place it in a beam and measure its performance, i.e., experimentally validate our 
ability to simulate precisely the passage of muons confined within a periodic lattice as 
they pass through energy absorbers and rf cavities. 

The experience gained from this experimental demonstration will provide important input 
to the final design of a real cooling channel. The successful operation of a section of a muon 
cooling cha nnel has been identified (most recently by the U.S. Muon Technical Advisory 



Committee 



157| ) as a key step in demonstrating the feasibility of a Neutrino Factory or 



Muon Collider. 
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B. Principle of the experiment 



Fundamentally, in a muon cooling experiment one needs to measure, before and afte r 



the cooling channel, the phase space distribution of a muon beam in six dimensions [15 
Such a measurement must include the incoming and outgoing beam intensities and must 
avoid biases due to the decay of muons into electron s wit hin the channel and due to possible 



contamination of the incoming beam by non- muons [159J. Two techniques have been consid- 
ered: i) the multi-particle method, in which emittance and number of particles in any given 
volume of phase space are determined from the global properties of a bunch; and ii) the 
single-particle method, in which the properties of each particle are measured and a "virtual 
bunch" formed off-line. The full determination of the covariance matrix in six dimensions 
is a delicate task in a multi-particle experiment, and the desired diagnostics would have to 
be developed specifically for this purpose; moreover, a high-intensity muon beam bunched 
at an appropriate frequency would need to be designed and built. For these reasons, the 
single particle method is preferred. The single-particle approach, typical of particle-physics 
experiments, is one for which experimental methods already exist and suitable beams are 
already available. 

In the particle-by-particle approach, the properties of each particle are measured in mag- 
netic spectrometers before and after the cooling channel fFigure EUj). Each spectrometer 
measures, at given z positions, the coordinates x, y of every incident particle, as well as 
the time. Momentum and angles are reconstructed by using more than one plane of mea- 
surement. For the experimental errors not to affect the measurement of the emittance by 
a significant factor, the rms resolution of the measurements must be s malle r than typically 
l/10th of the rms equilibrium beam size in each of the six dimensions 



C. Conceptual design 

Figure Wf\ shows the layout under consideration for MICE, which is based on two cells of 
the Feasibility Study II "Lattice 1" cooling channel. The incoming muon beam encounters 
first a beam preparation section, where the appropriate input emittance is generated by a 
pair of high-Z (lead) absorbers. In addition, a precise time measurement is performed and 
the incident particles are identified as muons. There follows a first measurement section, in 
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FIG. 46: (Color) Conceptual layout of MICE upstream spectrometer: following an initial time-of- 
flight (TOF) measurement, muons are tracked using detector planes located within a solenoidal 
magnetic field. Although in principle three x, y measurements as shown suffice to determine the pa- 
rameters of each muon's helical trajectory, in practice additional measurement redundancy will be 
employed; for example, a fourth measurement plane can be used to eliminate very-low-momentum 
muons that would execute multiple cycles of helical motion. A similar spectrometer (but with the 
time-of- flight measurement at the end) will be used downstream of the cooling apparatus. 

which the momenta, positions, and angles of the incoming particles are measured by means 
of tracking devices located within a uniform-field solenoid. Then comes the cooling section 
itself, with hydrogen absorbers and 201 MHz rf cavities, the lattice optics being provided by 
a series of superconducting coils; the pairs of coils surrounding each absorber have opposite 
magnetic fields ("bucking" solenoids), providing tight focusing. The momenta, positions, 
and angles of the outgoing particles are measured within a second solenoid, equipped with a 
tracking system identical to the first one. Finally, another time-of-flight (TOF) measurement 
is performed together with particle identification to eliminate those muons that have decayed 
within the apparatus. 

D. Performance 

Simulations of MICE have been carried out for a configuration including four tracking 
stations per spectrometer, each station consisting of three crossed planes of 500-micron-thick 
square-cross-section scintillating fibers (Figure l4"H|) . immersed in a 5T solenoidal field. Time 
of flight is assumed to be measured to 70 ps rms. As shown in Figure l4"H| measurement reso- 
lution and multiple scattering of the muons in the detector material introduce a correctable 
bias in the measured emittance ratio of only 1%. (For this study the effect of the cooling 
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FIG. 47: (Color) Schematic layout of MICE apparatus. 





FIG. 48: (Color)A possible MICE tracking-detector configuration. 

apparatus was "turned off" so as to isolate the effect of the spectrometers.) 

Figure loTH illustrates the muon-cooling performance of the proposed MICE cooling appa- 
ratus. The normalized transverse emittance of the incoming muon beam is reduced by about 
8%. The longitudinal emittance increases by about the same amount, thus the net cooling 
in six dimensions is also about 8%. These are large enough effects to be straightforwardly 
measured by the proposed spectrometers. 
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FIG. 49: (Color) Distribution of ratios of output to input six-dimensional emittance for 1000 sim- 
ulated experiments, each with 1000 accepted muons. The top figure shows the distribution of 
this ratio for the emittances as generated by simulation; the bottom figure, as "measured" in the 
simulated experiments. The curves are Gaussian fits to the points. 
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FIG. 50: (Color)Results from ICOOL simulation of MICE: normalized transverse (left) and longi- 
tudinal (right) emittances vs. distance. 



The CERN Neutrino Factory Working Group has studied a variant of the proposed MICE 
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cooling apparatus, in which 88-MHz rf cavities are employed in pl ace o f the 201-MHz devices 
(the 88- and 201-MHz designs have similar cooling performance) 16l| . Figure I5T1 ffrom the 
CERN study) elucidates further experimental issues. As shown in Figure l5Tk . for input 
emittance above the equilibrium emittance of the channel (here about 3500 mm-mrad), the 
beam is cooled, while for input emittance below equilibrium it is heated (and, of course, 
for an input beam at the equilibrium emittance, the output emittance equals the input 
emittance). Figure 15Tb illustrates the acceptance cutoff of the cooling-channel lattice; for 
input emittance above 6000 mm-mrad, the transmission probability falls below 100% due to 
scraping of the beam. Figure 15Tb shows the effect o f vary ing the beam momentum: cooling 
performance improves as the momentum is lowered 16J], as quantified here in terms of the 
fractional increase in the number of muons within the phase-space volume accepted by a 
hypothetical acceleration section downstream of the cooling channel. The goal of MICE 
includes verification of these effects in detail in order to show that the performance of the 
cooling apparatus is well understood. Subsequent running could include tests of additional 
transverse cooling cells, alternative designs, or emittance exchange cells. 

One critical aspect of this experiment is operation in the presence of backgrounds due 
to dark currents from the rf cavities. While it is possible to operate the experiment using 
comparatively low rf gradients, it would be highly desirable to produce cavities which would 
yield less dark current at higher gradients. This would permit more efficient use of the rf 
cavities and power supplies. We are trying to develop cavities with low dark currents. 



VII. CONCLUSIONS 



In summary, the Muon Collaboration is developing the knowledge and ability to create, 
manipulate, and accelerate muon beams. Our R&D program will position the HEP com- 
munity such that, when it requires a Neutrino Factory or a Muon Collider, we shall be in 
a position to provide it. A staged plan for the deployment of a Neutrino Factory has been 
developed that provides an active neutrino and muon physics program at each stage. The 
requisite R&D program is diversified over laboratories and universities and has international 
participation. 

The very fortuitous situation of having intermediate steps along this path, that offer a 
powerful and exciting physics program in their own right, presents an ideal scientific oppor- 
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FIG. 51: (Color)Simulation results for 88-MHz variant of MICE apparatus: a) output emittance vs. 
input emittance, with 45° line (dashes) superimposed; b) beam transmission vs. input emittance; 
c) cooling performance (see text) vs. input emittance for various beam kinetic energies (top to 
bottom: 140, 170, 200, 230 MeV). 

tunity, and it is hoped that the particle physics community will be able to take advantage 
of it. 
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